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Abstract. Let X be a homogeneous polynomial vector field of degree 2 on 
S'^. We show that if X has at least a non-hyperbolic singularity, then it has 
no limit cycles. We give necessary and sufficient conditions for determining if a 
singularity of X on S'^ is a center and we characterize the global phase portrait 
of X modulo limit cycles. We also study the Hopf bifurcation of X and we 
reduce the 16*'' Hilbert's problem restricted to this class of polynomial vector 
fields to the study of two particular families. Moreover, we present two criteria 
for studying the nonexistence of periodic orbits for homogeneous polynomial 
vector fields on of degree n. 



1. Introduction and statement of the main results 
A polynomial vector field X in is a vector field of the form 

d d d 

(1) X = p{x, y^^)-Q^ + Qi^^ y,z)-^ + R{x, y> ^)^' 

where P, Q and R are polynomials in the variables x, y and z with real coefficients. 
We denote n = maxjdegP, deg Q, deg i?} the degree of the polynomial vector field 
X. In what follows X will denote the above polynomial vector field. 

Let §2 be the 2-dimensional sphere {{x,y,z) g : x^ + y^ + z^ = 1}. A 
polynomial vector field X on is a polynomial vector field in such that restricted 
to the sphere defines a vector field on S^; i.e. it must satisfy the equality 

(2) xP{x, y, z) + yQ{x, y, z) + zR{x, y, z) ^ 0, 

for all points (x, y, z) of the sphere S^. 

Let / g R[a;,?/, 2:], where R[a;,y,z] denotes the ring of all polynomials in the 
variables x, y and z with real coefficients. The algebraic surface / — is an 
invariant algebraic surface of the polynomial vector field X if for some polynomial 

df df df 

K e M.[x,y, z] we have Xf = P-^ + Q-^ + R^ = Kf. The polynomial K is 

ox ay oz 

called the cof actor of the invariant algebraic surface / = 0. We note that since the 
polynomial system has degree n, then any cofactor has at most degree n — 1. 

The algebraic surface f — defines an invariant algebraic curve {/ = 0} n 
of the polynomial vector field X on the sphere if 
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points (x, y, z) G {/ = 0} H §^ we have that {x, y,z) A ( — — , — , 1 7^ 0, 



df df df 

(i) for some polynomial K e M.\x, y, z] we have Xf ~ P— — h Q— — h R-^ 

ax ay oz 

Kf, on all the points {x,y,z) of the sphere and 

(ii) the intersection of the two surfaces / = and is transversal; i.e. for all 

'a/ 9/ df 

^dx^ dy' dz ^ 
where A denotes the vector cross product in R'^. 

Again the polynomial K is called the cofactor of the invariant algebraic curve 

{/ = 0}ns2. 

Note that, if a curve {/ = 0} nS^ satisfies the above definition, then it is formed 
by trajectories of the vector field X. This justifies to call {/ = 0} n §^ an invariant 
algebraic curve, since in this case it is invariant under flow defined by X on S^. 

The equation of a plane in R'^ is given for ax + by + cz + d ~ 0. Any circle 
on the sphere lies in a plane ax + by + cz + d = 0, where we can assume that 

+ 6^ + = 1 and < —d < 1. If the invariant algebraic curve {/ = 0} fl 
is contained in some plane, then we say that {/ = 0} n §^ is an invariant circle of 
the polynomial vector field X on the sphere Moreover, if the plane contains the 
origin, then {/ = 0} H is an invariant great circle. 

Let U be an open subset of K'^. Here a nonconstant analytic function H : U ^M. 
is called a first integral of the system on U if it is constant on all solutions curves 
{x{t),y{t), z{t)) of the vector field X on U; i.e. H{x{t),y{t), z{t)) = constant for 
all values of t for which the solution {x{t),y{t), z{t)) is defined in U. Clearly H is 
a first integral of the vector field X on U ii and only if XH = on J7. If X is a 
vector field on S^, the definition of first integral on §^ is the same substituting U 
by i7n§2. 

In what follows we say that two phase portraits of the vector fields Xi and X2 on 
§^ are {topologically ) equivalent, if there exists a homeomorphism h : S'^ such 

that h applies orbits of Xi into orbits of X2 , preserving or reversing the orientation 
of all orbits. Similarly we define (topological) equivalence in the Poincare disc, see 
Section [5] for a definition. 

In [TT] we studied homogeneous polynomial vector fields of degree 2 on and 
we determined the maximum number of invariant circles when it has finitely many 
invariant circles. Moreover, we characterized the global phase portraits of these 
vector fields when they have finitely many invariant circles. Camacho in 1981 
proved some properties of this kind of vector fields, see also [10] . One of the results 
that can be found in [TT] and that will be necessary in this paper is given below. 

Proposition 1. Let X be a homogeneous polynomial vector field of degree 2 on M.^ . 
Then X is a polynomial vector field on §^ if and only if the system associated to it 
can be written as 

X = P{x,y, z) — aixy + a2y^ + asz^ + OiXZ + a^yz, 
(3) y = Q{x,y,z) ^ ~aix^ - a2xy + oqz^ + arxz + agyz, 

z = R{x, y, z) = —a/^x^ — a^y"^ — {a^ + ai)xy — a^xz — a^yz. 

The main results of this paper are the following ones. 

The next theorem characterizes the centers of the quadratic homogeneous poly- 
nomial vector fields on 



Theorem 2. Let X be a homogeneous polynomial vector field on of degree 2 
and let p — (0, 0, —1) be a singularity of X ; i.e. the system associated to X can be 
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written in the form ([3]) with — = 0. Then p is a center of X if and only if 
as — —04, a| + 0705 < and 04(0^ — a|) + 0102(05 + 07) = 0. 

Theorem [2] will be proved in Section ID 

In the next two results we provide sufficient conditions for the non-existence 
of periodic orbits (and in particular of limit cycles) for homogeneous polynomial 
vector fields on S'^ of arbitrary degree. For more details on limit cycles see [13] . 

We use the notation F{u, v) — F{2u, 2v, + — 1), where F e M.[x, y, z]. 

Theorem 3. Let X = {P,Q,R) he a homogeneous polynomial vector field on 
Then X has no periodic orbits if at least one of following conditions is satisfied. 

(a) The polynomial R does not change sign on and X does not have a 
periodic orbit passing through the point (0,0, 1). 

(b) The function ((P, Q), Vi?) |^_q does not change sign and X does not have 
a periodic orbit passing through the point (0, 0, 1). 

Now, we use the notation /C = aP + bQ + cR, Fi {u,v) — F{u + a,v + b,c— [an + 
bv)/c), F2{u, v) — F(a — {hu + cv)/a, u + b,v + c) and ^3(11, v) — F{u + a,b— [au + 
cv)/b,v + c) and F £ M.[x,y, z]. 

Theorem 4. Suppose that the homogeneous polynomial vector field X = (P, Q, R) 
does not have periodic orbits on §^ intersecting the great circle C determined by 
plane ax A- by + cz = with a^ + b^ + c^ — 1 . Then it has no periodic orbits if at 
least one of following conditions is satisfied. 

(a) The polynomial ICi does not change sign on M^. 

(b) The function ((Pi, Qi), V/Ci) |;ci=o ^'-'^^ ™^ change sign. 

(c) The polynomial IC2 does not change sign on . 

(d) The function (((52, P2), V/C2) Ij^^^o does not change sign. 

(e) The polynomial /C3 does not change sign on R^. 

(f) The function ((P3, P3), V/C3) lx;3=o ^"^'^ ™^ change sign. 
Theorems [3] and |4] will be proved in Section [S] 

The following theorem characterizes the phase portraits of quadratic homoge- 
neous polynomial vector fields on having at least one non-hyperbolic singularity. 

Theorem 5. Let X be a homogeneous polynomial vector field on §^ of degree 2. 

(a) // X has a linearly zero singularity, then its phase portrait is equivalent to 
the phase portrait of Figure QJ 

(b) If X has a nilpotent singularity, then its phase portrait in the Poincare disc 
is equivalent to one of the phase portraits of Figures\^ or0 

(c) // X has a center, then its phase portrait in the Poincare disc is equivalent 
to one of the phase portraits of Figures\B[ or\^ 

(d) IfX has a semi-hyperbolic singularity, then its phase portrait in the Poincare 
disc is equivalent to one of the phase portraits of Figures^ or[3 

We note that the unique non-hyperbolic singularity which does not appear ex- 
plicitly in the statement of Theorem[5]is a weak focus. Such singularities can appear 
in the quadratic homogeneous polynomial vector fields on §^ and such vector fields 
are studied in Theorems [8] and [9l Hence, from Theorem [5] it follows immediately. 

Corollary 6. Let X be a homogeneous polynomial vector field on §^ of degree 2. If 
X has at least a non-hyperbolic singularity on different from a weak focus, then 
X has no limit cycles on §^ . 
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Theorem [5] and Corollary [6] will be proved in Section [6l 

In the next theorem we give an upper bound for the number of singularities of 
homogeneous polynomial vector fields on of arbitrary degree. 

Theorem 7. Let X = {P, Q, R) be a homogeneous polynomial vector field on S'^ of 
degree n. If X has finitely many singularities on then X has at most 2(n^ — n+1) 
singularities on . 

Theorem [7] will be proved in Section [71 

The following theorem characterize modulo limit cycles the phase portraits of 
quadratic homogeneous polynomial vector fields on S'^ having non-degenerate sin- 
gularities. 

Theorem 8. Let X be a homogenous polynomial vector field on §^ of degree 2. // 
all singularities of X are non-degenerate, then its phase portrait in the Poincare 
disc is topologically equivalent modulo limit cycles to one of the phase portraits 
giving in Figures fTSl 

Theorem |8] will be proved in Section [H 

The next result characterizes the families of quadratic homogeneous polynomial 
vector fields on §^ having a Hopf bifurcation. 

Theorem 9. Let X be a family of homogeneous polynomial vector fields on 
of degree 2 having a Hopf bifurcation. Then, doing a orthogonal linear change of 
variables, the system associated to X can be written as 

X = aixy + 02?/^ + asj/z, 

ij = —aix^ — a2xy + ajxz + a^yz, 

z = -{ar, + ar)xy - asy^, 

with fli 0, a2(a5 + a-j) — aiflg ^ and —a^ar < 0. 
Theorem ini will be proved in Section O 

We have reduced the possible existence of limit cycles to two families; more 
precisely, the families ((38)) and ([40]) . In Theorem [9] we have proved that family ([40]) 
has for some values of the parameters limit cycles. Finally in the last section we 
conjecture that family ([38]) has no limit cycles. 

2. Poincare disc 

Let X be a homogeneous polynomial vector field on S"^ , then the differential sys- 
tem associated to it is invariant with respect to the change of variables (x, y, z, t) 
(— —y, —z, —t) if its degree is even, or with respect to {x, y, z, t) i— > {—x, —y, —z, t) 
if its degree is odd. Thus, in particular the phase portrait of X at the northern 
hemisphere of is symmetric with respect to the origin of the sphere to the phase 
portrait at the southern hemisphere with the time reverse if the degree of X is 
even, or with the same time if the degree of X is odd. We now project the northern 
hemisphere of orthogonally onto the plane 11 containing the equator of i.e. 

The orbits of X on the northern hemisphere of are mapped onto certain 
curves of the unit disc of 11. We call this unit disc, together with the corresponding 
induced phase portrait, the Poincare disc. 

Now we consider the homogeneous polynomial vector field X = {P,Q,R) of 
degree m. We identify as the tangent plane to the sphere at the point 
p — (a, 6, c), i.e with the plane ax + by + cz — 1 = 0, where + b'^ + c^ — 1. 
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Suppose that c 7^ 0, then we denote the points of as {u + a,v + b,c — {au + 
bv)/c). Let TT : ^ n {ax + by + cz > 0} be the difFeomorphism given 
by 7r(M, v) ~ \c\/\/X {x ^ u + a, y ^ v + b^ z — c ~ (au + bv)/c), where X — c^{l + 
+ w^) + (au + bv)'^. That is, tt is the inverse map of the central projection 
TT-^ ; §2 n {ax + by + cz>0}^M.'^ defined by 

(4) n-\x,y,z)^ 



ax + by + cz ax + by + cz ax + by + cz ^ 

The homogeneous polynomial system X on §^ becomes, through the central pro- 
jection TT^^, the differential system 

u = ^{P-{u + a)(aP + bQ + cR)), 

r 

V = i^{Q~[v + b){aP + bQ + cK)), 
|c| 

on the plane M? . Here F — F {tt{u,v)). Since X is a homogeneous polynomial 
vector field of degree m we have that 

K 
Vx 

K 
Vx 



P - {u + a){aP + bQ + cR) 
(^Q-{v + b){aP + bQ + cR)'^ 



where F = F{u + a, 1; + 6, c — [au + bv)/c). If t denotes the independent variable 
in the above differential system, then this system becomes polynomial introducing 
the new independent variable s through ds — {^/X/\c\)^~"^dt, i.e. 

ii = V{u,v) = P-{u + a){aP + bQ + cR), 
V ^ Q(u,v)^Q-{v + b){aP + bQ + cR). 

Now in the case a ^ the central projection ^ induces a polynomial vector 
field on the plane ax + by + cz — I ~ determined by 



(6) 



it = ■P{u,v) = Q - {u + b){aP + bQ + cR), 
V = Q{iL,v) = R-{v + c){aP + bQ + cR), 



where F — F{a — (bu + cv)/a, u + b,v + c); and in the case 5 7^ we obtain 
(7) 



u = Viu^v) = P - {u + a){aP + bQ + cR), 
V ^ Q{u,v)^ R- {v + c){aP + bQ + cR), 



where F — F{u + a,b— {au + cv)/b, v + c). 

The proof of the next two propositions can be found in I 17| . 

Proposition 10. Let X ~ (P, Q, R) be a homogeneous polynomial vector field of 
degree m on . Then the planar vector field induced from X through the central 
projection ([4]) on the plane ax + by + cz ~ 1 — with a^ +b^ + c^ = 1 has degree m 
if and only if the great circle {ax + by + cz ~ 0} CiEi^ is an invariant circle of X . 

In general we will consider the case (a, b, c) ~ (0, 0, —1). Thus ([5]) becomes 

.„s ii = V{u,v) = P{u,v,—1) + uR{u,v,—l), 

V = Q{u,v) = Q{u,v,-l) + vR{u,v,~l). 
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In particular if X is a homogeneous polynomial vector field on §^ of degree 2 
associated to ([3]), then ^ becomes 

u — — a^u — flsw + a^u^ + (ai + aG)uv + aiv^ 
/„s —a^v? — {a^^a'j)u^v~a^uv^, 

V — flg — ayM — agw — aiu^ + (03 — a2)uw + agu^ — 
a^u'^v — (as + ai)uv'^ — agv^ . 

Proposition 11. Let X be a homogeneous polynomial vector field on §^ of de- 
gree m, and let f G y, z] be a homogeneous polynomial of degree n such that 
f{x,y,0) ^ 0. Then f = is an invariant algebraic surface of X if and only if 
f(u^ V, —1) = is an invariant algebraic curve of the polynomial planar vector field 
given by ([8]) induced from X through the central projection (j4|) . 

The next results jointly with their respective proofs can be found in Camacho 
[2] or in [HI- 

Proposition 12. Let X be the vector field associated to and let si, — Si G §^ 
be saddle points of X with a common separatrix I, then I is contained in a great 
circle passing through Si and —Si. 

3. Stereographic projection 

We identify as the plane ax + by + cz = with a^ + b^ + = 1. Suppose 
that c 7^ 0, then we denote the points of K.^ as {u,v, —{an + bv)/c). Let tt : R-^ ^ 
§^\{(a, 6, c)} be the diffeomorphism given by 7r(u, w) — 1/X{x — aX—2c'^{a—u), y = 
bX-2c'^{b-v), z = cX-2c{c^+au+bv)), where A = c^{l+u'^+v'^) + {au+bv)'^ . That 
is, TT is the inverse map of the stereographic projection tt^^ : \ {(a,5, c)} —^ 
defined by 7r~^(a;, z) = 



(10) 



X — a{ax + by + cz) y — b{ax + by + cz) z — c{ax + by + cz) 



1 — {ax + by + cz) 1 — (ax + by + cz) 1 — [ax + by + cz) 

Through the stereographic projection tt^^ the homogeneous polynomial system X 
on §^ becomes the differential system 

u = :^{P+{u~a){aP + bQ + cR)), 
V = —{Q+{v~b){aP + bQ + cRj), 

on the plane M^. Here F — F {Tr{u,v)). Since X is a homogeneous polynomial 
vector field of degree m we have that 

' P+iu-a){aP + bQ + cR) 



2c^>^ 



-{Q + {v- b){aP + bQ + c^)) 



where F = F{aX - 2c^{a -u),bX- 2c^{b -v),cX- 2c{c^ +au + bv)). If t denotes 
the independent variable in the above differential system, then this system becomes 
polynomial introducing the new independent variable s through ds = A^^™/ (2c^)dt, 
i.e. 

. s ii = V{u,v) = P+{u-a){aP + bQ + cR), 

^ ^ V = Q{u,v) = Q + {v-b){aP + bQ + cR). 
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Now the dot denotes derivative with respect to the variable s. In the cases a ^ 
the stereographic projection (llOp induces a polynomial vector field on the plane 
ax + by + cz — determined by 

, . u = V{u,v) = Q + {u-b){aP + bQ + cR), 

^ ' V = Q{u,v) = R+{v-c){aP + bQ + cR), 

where F = F{a\i — 2a{a^ + bu + cv),b\i — 2a?{b — u),c\i — 2a?{c — w)), with 
Ai = a^(l + V? + v^) + {bii + cti)^; and in the case & ^ we get 

,-^„x ii = 'P{u,v) = P + {u- a){aP + bQ + cR), 

^ ' V ^ Q{u,v) = R+iv-c){aP + bQ + cR), 

where F = F{a\2 - 2b'^{a - u),b\2 - 2b{b'^ + au + cv),cX2 - 2b'^{c - v)), with 
A2 = 6^(1 + + v^) + {au + cv)^. 

The planar vector field induced by the stereographic projection (fTO|) on the plane 
ax + by + cz = with + 6^ + = 1 will be called the stereographic projection of 
the vector field X = (P, Q, R) o,t the point (a, b, c) € 

In general we will consider the case (a, 6, c) = (0, 0, 1). Thus ((TT|) becomes 

, , u = P{2u,2v,u^ + v'^ ~l) + uR{2u,2v,u'^ + v'^ -1), 

^ ' V = Q(2u,2v,u'^ + -l)+vR{2u,2v,u'^ + -1). 

The proof of the next proposition can be find in [12| or in |17| . 



Proposition 13. Let X = (P, Q, R) be a homogeneous polynomial vector field of 
degree m onS^ . Then the planar vector field induced by the stereographic projection 
pop on the plane ax + by + cz — with a^ + b^ + c^ — 1 has degree 2m. Moreover, 
it has degree 2m — I if {a, b, c) ~ (0, 0, 1) is a singularity of X on S^. 

4. Classification of centers of the quadratic homogeneous 
polynomial vector fields on §^ 

The classification of the centers of the quadratic homogeneous polynomial vector 
fields on §^ is equivalent to the classification of the centers of system ^ induced 
from ^ by the central projection fl]). 

Proposition 14. Let X be a homogeneous polynomial vector fields on §^ of degree 
2. Then the system associated to X can be written in the form ([3]) with a^ = ag = 0. 

Proof. By the Poincare-Hopf theorem (see [U]), a system on §^ always have at least 
one singularity. Therefore, we can suppose that ^ has a singularity at (0,0, —1), 
because we can do a rotation of SO{3) which preserves all the properties of X. 
Hence, without loss of generality (0, 0) is a singularity of This means that we 
can suppose that in ([3]) 03 = ae = 0. □ 

By Proposition [Til ^ becomes 

ii~P{u,v) — ~a4U — a^v + aiuv + a2v'^ — a4U^— 

{a^ + ar)u'^v - asuv"^, 
v — Q{u,v) = — ayu — flgw — aiu^ — a2uw — a4M^w— 

(05 + ai)uv'^ — asv^. 



Let X — (P, Q). Consider the linear part of X on (0, 0), i.e. 
DX{0,0) = 



-04 —05 

-07 -as 
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Figure 1. The phase portrait of system ([3]) for as = 04 = as = 
ag = 07 = ag = 0. 

and denote by tr{DX{Q,0)) = —04 — ag and det{DX{0,0)) — a^ag — a^aj its 
trace and its determinant. If (0,0) is a center of (fTSl) . then after a linear change 
of variables and a rescaling of the time variable, it can be written in one of the 
following three forms: 

r ^ -s + P2{r,s), s = r + Q2{r,s). 
called a linear type center, 

r^s + P2{r,s), s^Q2{r,s). 

called a nilpotent center, 

f^P2{r,s), s^Q2{r,s). 

called a linearly zero center, where P2, Q2 are polynomials with terms bigger or 
equal than 2. 

We say that p = (0,0,-1) is a linear type center oi X if (0,0) is a linear type 
center of (fT5|) . In the same way we define nilpotent centers and linearly zero centers 
of X. 

Proposition 15. Let X he the vector field associated to system ()15|) . //a4 = 05 = 
07 = ag = 0, i.e. DX{{),Q) is the zero matrix, then (0,0) is not a center of X . 
Moreover, the phase portrait of ^ is equivalent to Figure[^ 

Proof. By hypothesis system ([3]) becomes 

X = y{aix + a2y), y = -x{aix + a2y), i = 0. 

Therefore, the phase portrait of ([3]) is equivalent to Figure [T] Hence (0,0) is not a 
center oi X. □ 

Proposition 16. The origin of (|15p cannot be a nilpotent center, i.e. if 
tr(DX(0,0)) = det(L»X(0,0)) = and \a4\ + jasj + |a7| + jagl ^ 0, then (0,0) 
is not a center. 

Proof As tT{DX{0,0)) = dct(DX(0,0)) = 0, i.e. ag = -a4 and aj + 0705 = 0, we 
suppose first that a4a7 7^ 0. Therefore, as = — a4/a7. Doing the change of variables 
It = r — s/a4, V = ajs/a^ system (jl5|) can be written as 

(16) r^s + P2ir,s), s = g2(r,s). 

Let (p{r) = air + a2r'^ + ■ ■ ■ the solution of s + P2{r, s) = 0. Therefore, ai = 
and a2 = {aia^ + 0207) /aj. If a2 — 0, then in the above system (0,0) cannot 
be a center because in this case it is not an isolated singularity. Now, if 0:2 7^ 0, 
substituting (p in ip{r) = Q2{r,(p{r)) we obtain that ip{r) = /32?'^ + where 
P2 = {{0,4 + a^){a4ai + aja2))/ {a4a7). Hence, as 02 7^ 0, it follows, by Theorem 67 
of page 362 of [1,, that (0,0) is not a center. 
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Now, we have to study the case al + asoy = 0, 0407 = and |a4| + jasj + layj 7^ 0. 
Therefore, we have two possibiHties 04 = 05 = 0, ay ^ and — — 0, ^ 
for system ([T5|) . But, as these two possibihties are equivalent by the change of 
variables {u,v) 1-^ {v,u), we need to study only the last one. Therefore, after a 
rescaling of time given by dr = —a^dt in (|15p . we obtain a system of the form 

(17) U = S + P2(u,w), V^Q2{u,v). 

Note that ai 7^ in (fTT)) . otherwise (0,0) is not an isolated singularity. As in 
the previous case, we have that (p{u) = is the solution of u + P2{u,v) — and 
ip{u) = Q2{u, f{u)) — {aiu )/a5 + • • • . Hence, by the same argument (0, 0) is not 
a center. □ 

Before proving Theorem [5] we do one remak. 

Remark 17. Let U be an open subset o/M^, F : U R be an analytic function 
which is not identically zero on U and X be a polynomial vector field associated to 
a system of the form 

(18) x^P{x,y), y^Q{x,y), 

where P and Q are polynomials in the variables x and y with real coefficients. The 
function F is an integrating factor of this polynomial system on U if one of the fol- 

d(FP) d(FQ) 

lowing three eguivalent conditions holds on U — — , div{FP, FQ) ~ 

ox ox 

0, XF — —F div(P, Q). As usual the divergence of the vector field X is defined by 
div(X)=div(P,Q) = g + ^. 

The first integral H associated to the integrating factor F is given by 

(19) H{x, y)^ J Fix, y)P{x, y)dy + h{x), 

where h(x) is chosen in order that it satisfies — — = —FQ. Note that — — = FP, 

ox oy 
so that XH = 0. The function H is single-valued, if U is simply connected and 
F^O. 

Conversely, given a first integral H of the system associated to X we always can 

OH dH 
find an integrating factor F such that — — = FP and — — = —FQ. 

Oy ox 

Let V : U R be an analytic function which is not identically zero on U . The 
function V is an inverse integrating factor of the polynomial system il8\) on U if 

^ ^ dx ^ ^ dy ^dx ^ dy^ ' 

We note that {V — 0} is formed by orbits of system US]} and F = 1/V defines on 
U \ {V — 0} an integrating factor of I118\) . 

Proof of Theorem[B We have that if tr{DX{0,0)) = and det{DX{0,0)) > 0, 

1. e. flg = —04 and al + 0705 < 0, then the origin can be a linear type center or a 
weak focus. Note that in this case 0507 7^ 0. Hence, doing the change of variables 

u = —\J—{a\ + ara^)r/aj + Ois/ai, v = s and introducing the new time variable 

T through dr = (a| + 0705) dt system (|15p can be written as 

(21) r = -s + P2{r,s)+P3{r,s), s = r + Q2{r, s) + Q^ir, s), 
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where P2, P3, Q2 and Q3 are respectively the homogeneous parts of degree 2 and 
3 of the above system. We have that the origin of this system is a center if and 
only if its Lyapunov coefficients Vt are zero. Therefore, we have to determine the 
Lyapunov coefficients Vi that can be obtained from 

(22) ^|^+o§-f:^'(-'+^'r^^ 

1=1 

where H{r, s) — l/2(r^ + s^) + Hj{r, s) and Hj are homogeneous polynomial of 

3=3 

degree j . 

We will show that if Vi — 0, then Vk — for all k E N. From the homogeneous 
part of degree two, three and four of (P^ . we obtain the following system 
(23) 

dH2 , dH2 
hr— V P2^ yQ2^^ = 0, 



dm ^ dHi dir^ dll3 dH2 0U2 . 

-s^ hr— — +P2-^ +Q2-^ +^3^ ^^3-^— = V^ir-^ + s). 

or OS Or OS or Os 

Solving system we obtain that Vi — ((05 — a^){a4^{a\ — a^) + 0102(05 + 

07)))/(8o7(— al — osfly)-^/^). Now, as a\ + osar < 0, we have that 05 7^ 07. There- 
fore, Pi = if and only if a4^{a\ — a^) + 0102(05 + ay) = 0. Hence, we distinguish 
the following four cases. 

Case 1: oi = 04 = and 02 ^ 0. In this case system ()15p becomes 

(24) u = —a^v + 02W^ — (05 + ai)u^v, v — —a-ju — a2UV — (05 + a-;)uv^ . 

This system is invariant with respect to the change of variables {u, v, t) 1— > (— u, w, — t), 
i.e. it is symmetric with respect to the straight line m = 0. Therefore, in this case 
(0, 0) is a center. 

Case 2: 02 = and oi 7^ 0. Then 04 = 0. In this case system (fT5|) becomes 

(25) li = —0517 + aiuv — (05 + a'j)u'^v, v — —a'ju — aiu^ — (05 + a'j)uv'^ . 

We have that system ([25l) becomes (|24|) by the change of variables (m, z;) i-^ {v, u). 

Therefore, this case is equivalent to Case 1. 

Case 3: 01=02 = 0. In this case system ([T5|) becomes 

, - u — — 04?/ — 051; — 04?/'^ — (05 + a7)w^u + 04uu^, 

i) = —a^u + a4V — a4U^v—{a5+a7)uv'^+a4V^. 

This system has a first integral given by H{u, v) — ((07 + 05)1*^ — 2a4ww + a5)/ ((05 — 
07)^^ + 2a4uv + 2a5t;^ + 05). As H is defined in (0, 0), it follows that in this case 
(0, 0) is a center. 

Case 4- ai02 7^ 0. In this case Vi = implies that 05 = (04(02— o^)— 010207)7(0102). 
Note that — (04 + O5O7) > becomes ((oia4 + 0207)(ai07 — a204))/(oi02) > 0. Now, 
this case is equivalent to Case 1 doing the orthogonal linear change of variables 
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where the a^'s denote the coefficients of X in the variables {x,y,z) satisfying the 
hypothesis. In fact, denoting by ai the coefficients of system ([3|) in the variables 
(i, z), we have that di = 0,3 — 0,4 — ~ iig, ~ 0, 02 = —yui^+af, 05 ~ 
(aioy — a2a4)/ai, 07 = —(0104 + a2aj)/a2. Note that 02 7^ and 0507 < 0. This 
completes the proof of the theorem. □ 

5. Two CRITERIA FOR DETERMINING THE NONEXISTENCE OF LIMIT CYCLES FOR 
HOMOGENEOUS POLYNOMIAL VECTOR FIELDS ON 

Let X — {P, Q, R) be a homogeneous polynomial vector field on §^ of degree m. 
The next proposition will be necessary for proving the results of this section, see 
its proof in [11] . 

Proposition 18. Let f{x, y,z) = + + — 1 and let X he a homoge- 
neous polynomial vector field on Then f is a first integral of X in 'E? , i. 
e. P{x, y, z)x + Q(x, y, z)y + R{x, y, z)z = for all (x, y, z) S R'^. 

Let (a, 6, c) G S^, we can suppose that X does not have a limit cycle such 
that (a, 6, c) belongs it, because we can do a rotation of 50(3) which preserves all 
properties of X . Hence, the study of the existence of limit cycles of X on is 
equivalent to study the existence of limit cycles of the planar vector field induced by 
X through the stereographic projection pO)) at the point (a, 6, c). We will consider 
the cases when (a, 5, c) is the point (0, 0, 1), i.e. the planar vector field determined 
by system ([T4|) . We denote this planar vector field by 

X = (P, g) = (P + uR, Q + vR). 

Here F{u, v) ^ F{2u, 2v , u"^ + v'^ - 1) with F eR[x,y,z]. 

Now we state a result that can be found in [8J and that will be necessary for the 
proof of Theorem [3l 

Theorem 19. Let (V, Q) be a vector field defined in an open subset U of R^, 
{a{t), l3{t)) a periodic solution of {V,Q) of period t, K : U — > M. a map 
such that K{a{t), P{t))dt 7^ 0, and f = f{x,y) a solution of the linear 

partial differential equation + Qtt" = ^ f- Then the closed trajectory 7 — 

ox ay 

{{a{t) , I3(t)) G U : i e [0, r]} is contained inT, — {{x, y) £U : f{x, y) = 0}, and 7 
is not contained in a period annulus of {V, Q). Moreover, if the vector field (P, Q) 
is analytic, then 7 is a limit cycle. 

Proof of Theorem\^ Consider the function f(u, v) — + + Using the previous 
notation we have that 

Pw- + Qw- = 2mP + 2u^R + 2vQ + 2v^R = 2uP + 2vQ + 2{ip- + v^)R. 
ou ov 

Note that, by Proposition [TH 2uP + 2vQ + {u^ + v"^ - l)R = 0. Therefore, 

+ - -(w' + - l)R + 2{u^ + v^)R = R{u^ + v^ + l)^ Rf. 
ou ov 

Therefore, / = is an invariant algebraic curve of X and R is its cofactor. 

Using the notation of Theorem [T9l we assume that 7 is a periodic orbit of X. 
Now, if R does not change sign in R^, then J^Rdt ^ 0, and since {{u,v) S R^ : 

f{u,v) = 0} = 0, by Theorem [19] X does not have periodic orbits. Hence, X does 
not have periodic orbits on S'^. This proves statement (a). 
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Now we will prove statement (b) . Note that 

Therefore, if ((P,Q),Vi?) does not change sign, then the vector field X is 

transversal to the curve i? = 0, so there are no periodic orbits of X intersecting the 
curve R = 0. Hence, if X has a periodic orbit 7, then 7 C \ {B. = 0}. Thus, 
J^Rdt ^ 0, and so by Theorem [T9] we get 7 C G : f{u,v) = 0} = 0. 

Therefore, X does not have periodic orbits, i.e. X does not have periodic orbits on 
S\ □ 
From the proof of Theorem [3] we have: 

Corollary 20. The vector field X has the invariant algebraic curve f{u,v) = 
u'^ + + 1 = with cofactor Kf — R. 

We see in the beginning of this section that the study of the existence of limit 
cycles of homogeneous vector fields X on is equivalent to study the existence 
of limit cycles of the planar vector fields induced by X through the stereographic 
projection at the point (a, 6, c). But sometimes it is necessary to do a rotation in 
order that (a, 6, c) does not belong to a limit cycle of X. In fact, if (a, 6, c) belongs 
to the limit cycle 7 then, as the phase portrait of X is symmetric with respect to 
the origin, we have that —(a, 6, c) G —7 = {—{x, y, z) : [x, y, z) G 7}. If degree of 
X is even, by Proposition 4 of [H], 7 7^ —7. Therefore, if the degree of X is even, 
it is not necessary to do a rotation. Now if degree of X is odd, by Proposition 6 of 
[T^ we know that the case 7 = —7 is possible. 

Now, we will consider the the planar vector fields Xi, i = 1, 2, 3, determined by 
systems (O, dH) and ([42]), i.e. 

Xi = iPi,Qi) = {Pi-{u + a)ICi,Qi-iv + b)K:i), 

X2 = [P2.Q2) = {Q2-[u + b)jC2,R2-[v + c)iC2). 
X3 = {P3,Q3) = {P3-{u + a)IC3,R3-{v + c)IC3). 

Here K. = aP + bQ + cR, Fi{u, v) = F{u + a,v + b,c — {au + bv)/c), F2{u, v) — 
F{a — {bu + cv)/a, u + b,v + c) and F^iu, v) — F{u + a,b — (au + cv)/b, v + c) and 
FeR[x,y,z]. 

Proof of Theorem First we will prove statement (a) . Consider the function 

. , , , x9 / ,^9 ( au + bv^ 

fi{u,v) = {u + aY + {v + by + [c- 



c 

2ab P + 2 

— T^UV H V - 



Using the previous notation we have that 

— \-wi^r- — ^1-5 — t-Qi^ 
au ov ou ov 



Pi^+Qi^ = Pi^ + Qi^-^iCi{fi-l)-\iCi{au + bv) 



= -2^1 (/i - 1) + ^{c^uPi + c^vQi - c{au + bv)Ri). 

Note that, by Proposition [TH {u + a)Pi + {v + b)Qi + (c - {au + bv)/c)Ri = 0, i.e. 

c^uPi + c^vQi - c(au + bv)Ri = -c^JCi. Therefore, A^A + Oi^A = -2^i(/i - 

ou ov _ 

1) — 2/Ci = — 2/Ci/i. Hence, /i = is an invariant algebraic curve of Xi and — 2/Ci 
is its cofactor. 
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Using the notation of Theorem [19] we assume that 7 is a periodic orbit of Xi . 
Now, if K-i does not change sign in M^, then J^ICidt ^ 0, and since {{u,v) G : 

fi{u,v) = 0} = 0, by Theorem[T9]Xi does not have periodic orbits. Hence, X does 
not have periodic orbits on This proves statement (a). 
Now we will prove statement (b) . Note that 

Therefore, if ((Pi , Qi), V/Ci) | 

K.i=o '^ot change sign, then the vector field Xi 
is transversal to the curve /Ci = 0, so there are no periodic orbits of Xi intersecting 
the curve ACi = 0. Hence, if X has a periodic orbit 7, then 7 C ]R^\{ACi = 0}. Thus, 
J^ICidt ^ 0, and so by Theorem [H] we get 7 C {{u,v) G : fi{u,v) = 0} = 0. 

Therefore, Xi does not have periodic orbits, i.e. X does not have periodic orbits 
on §2. 

The proof of the other statements is similar replacing /i by 

/2 = {u + by + {v + cy + [a 

in the statements (c), (d) and by 

h[u,v) = [u + af + {v + cf + [b 



in the statements (e), (/). □ 
^From the proof of TheoremlU using the previous notation, we have the following 
corollaries: 

Corollary 21. The algebraic curve f i = is an invariant algebraic curve of vector 
field Xi with cof actor Kf. ~ — 2/Ci, for i — 1,2,3. 

Corollary 22. In the assumptions of Theorem^ if C is a periodic orbit, then C 
is the unique periodic orbit of X on S^. 

We say that an orbit 7 of a homogeneous polynomial vector field on is convex 
if for any p € 7 the great circle tangent to 7 at p intersects 7 only at p, or 7 is 
entirely contained in the great circle. 

Proposition 23. Let X — (P, Q, R) be a quadratic homogeneous polynomial vector 
field on §^ . Any great circle in §^ has at most four tangencies with the orbits of X 
or is formed by orbits of X . 

Proposition 24. Let X be a homogeneous polynomial vector field on §^ of degree 
2. Any periodic orbit of X on §^ is convex. 

Propositions [23] and [M] also can be found in Camacho [5] or [T7] . 

By Proposition [24] we have that in the case that X has degree 2 the hypotheses 
of Theorem [4] is not necessary, i.e. we can suppose that all periodic orbits of X are 
contained on the hemisphere determined by the plane ax -\- by + cz = such that 
(a, 6, c) belongs it. 
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6. Phase portraits for homogeneous polynomial vector fields on 
of degree 2 having some non-hyperbolic singularity 

Let p be a singularity of a homogeneous polynomial vector field X on of degree 
2. By Proposition [T31 we can suppose without loss of generality that p = (0, 0,-1), 
then the classification of singularities of X is equivalent to the classification of 
singularities of (fTS]) at the origin. 

Using the notation of the previous section we say that p is a hyperbolic singularity 
of X if the eigenvalues of DX{0, 0) have real part distinct of zero. We say that p is 
a non-degenerate singularity if det(DX(0, 0)) = 0403 — 0507 7^ 0, a semi-hyperbolic 
singularity if tr(_DX(0,0)) = —04 — ag 7^ and det(-DX(0, 0)) = 0, a nilpotent 
singularity if tr(i:>X(0, 0)) = det{pX{0,0)) = and |a4| + \ar,\ + \a7\ + lagj ^ 0, 
and linearly zero singularity if DX{0,0) is the null matrix. 

The next result will be necessary for proving the results of this section and its 
proof can be found in [llj . 

Proposition 25. Let X be a homogeneous polynomial vector field of degree 2 on 
Suppose that X has at least one invariant circle on then X has no limit 
cycles. 

Proof of Theorem\^ Statement (a) is exactly Proposition 1151 Therefore, we have 
to prove only statements (&), (c) and (d). 

In what follows, by Proposition [TH we can suppose that the system associated 
to X = (P, Q, R) is in the form ([3]) with as = ag = 0. Therefore, if i? ^ then the 
equator of i.e. S"'^ is not invariant by X, because z is not a factor of R. Hence, 
in this case, on the Poincare disc associated to X we represent §^ by a circle formed 
by dots. 

In the case that X has a nilpotent singularity, to determine its phase portrait, 
from the proof of Proposition [TBI it is sufficient to study the systems 



,2 



X = aixy + a2y + a4a;z yz, 

(27) y = = —aix — a2xy + a^xz — a/^yz, 



-a^x"^ + 04?/^ — I a? ^ I xy. 



a, 
ay 



with a4ay 7^ 0, and 

(28) X = aixy + a2y'^ + a5yz, y = -aix^ - a2xy, z ^ -a^xy, 

with as 7^ 0. Now system ([27| is equivalent to system (p8|) by the orthogonal linear 
change of variables 




/ - /V \ 



ay a4 







X 

V 



V 1 y 

where Oi are the coefficients of system (P7)) . Denoting by hi the coefficients of 
system (1271) in the variables (x,y, 5), we have that ai = (aia4 + a2ay) I \J a\-V ay, 
0,2 = {0.20,4 — aiay) / a1 + fly, 0,5 ~ —(04 + ay)/ay- Note that 7^ 0. Therefore, 
we have to study only system (j28p or equivalently system (jl7p induced from it by 
the central projection. We distinguish two cases. 
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Figure 2. Phase portrait of Case 1: aia2 7^ 0. 



Case 1: ai ^ 0. By the proof of Proposition [TBI we have that v = (p{u) = 
is the solution of u + P2{u,v) = 0. Hence, ip{u) = Q2{u,ip{u)) = oiu^/as and 

dP2 , , , dQl , , ^^ a2 



dv 



iu,ip{u)) = —u. 

as 



Then, by Theorem 67 of page 362 of ^ it 



does not have 



follows that (0, 0) is a cusp. 

In the case 02 7^ the vector field X determined by system 
singularities on S^. Moreover, X has four singularities (0, 0, ±1), (0, a^/{\Ja\ + a|), 
—a2/{yja\ + flg)) and (0, — a5/(v^af+lif ), a2/( v^of+lif)) corresponding to the 
singularities (0,0) and (0,05/02) of (fT7|) . We have that the eigenvalues associ- 
ated to (0,05/02) are (— oi ± yja\ — 4(o| + Og))/2o2. Thus, (0,05/02) can be a 
node or a focus of (fT7|) . Now, by the orthogonal linear change of variables 











0,2 





05 



05 



2 

02 




the system associated to X becomes 

0102 



(29) 



P{x,y,z) = 



V4 + 



0105 



- «5 

ai02 



=£2 



- 02^^ + 05x1, 



Note that in these coordinates the singularity (0, a^l{\J + Og), — 02/(\/o| + Og)) 
correspond to the point (0,0, —1). Hence, using the notation of SectionlHl we have 
that R{u,v) = R{2u,2v,v? + — 1) = Aaia^/ yja^ + a^v? . Thus, R does not 
change sign and by Theorem [31 X does not have periodic orbits on So the 
phase portrait of X on the Poincare disc is equivalent to one of Figure [H 

In the case 02 = the vector field X has four singularities given by (0, 0, ±1) 
and (0, ±1, 0). The eigenvalues associated to the last two singularities are 0, (oi ± 

\J ai — 4a5)/2 and 0, — (oi ± \J a\ — Aa\)/2, respectively. Therefore, these singu- 
larities can be nodes or foci. Now, substituting 02 = in the above change of 
coordinates we obtain that in this case the system associated to X is equivalent to 
system ([29)1 with 02 = 0. Therefore, by the previous argument, we have that X 
does not have limit cycles on S^. Thus, the phase portrait of X on the Poincare 
disc is equivalent to the one of Figure [H 

Case 2: oi = 0. In this case we have that (fT7|) has the straight line of singularities 
V = Q and, so X has an invariant circle formed by singularities determined by the 
intersection of the plane y = with Moreover, H{x, y, z) = a^y — a2Z is a first 
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Figure 3. Phase portrait of Case 1: ai ^ and a2 — 0- 




(«) {b) 
Figure 4. Phase portrait of Case 2. 

integral of X (note that 05 7^ 0) and X has other two singularities (0, a^, — 02)7(02 + 
a|) and (0, — 05, a2)/(a2 + 0,5)- Then its phase portrait is equivalent to the one of 
Figured] (a) if 02 7^ 0, or to the one Figured] (&) if 12 = 0. Hence, statement (b) is 
proved. 

Now, we prove statement (c). By Theorem |2l we can suppose that X is in the 
form ([3]) with = = 0, ag = —04, al+aja^ < and 04(0^ — a|)+aia2(a5+a7) = 
0. Moreover, by the proof of Theorem |2l for characterizing the phase portrait of 
X we have to study only systems ([^ and (pS)) . Therefore, we distinguish the 
following two cases. 

Case 3: ai — Q and 02 ^ 0. Then 04 = 0. This case correspond to system ([24| . 
First we suppose that aj 7^ —05. In this case X has only two singularities on 
given by (±1,0,0) having eigenvalues {0, — (02 ± ^^a^ — 4(0^ + a507))/2} and 
{0, (02 ± \/a^ — 4(ay + a5a7))/2}, respectively. The system induced by X through 
the central projection, i.e. system has a center in (0, 0) and another singularity 
(0,05/02) with eigenvalues ±-\/— (oj + a^ia^ + 0507)702. Note that 

(30) (05 + 0507) (07 + 0507) — (05 + 07)^0507 < 0, 

because 0507 < and 05 + 07 7^ 0. Therefore, if the singularities on §^ are sad- 
dles, then (0, 05/02) is a center because system is invariant with respect to the 
change of variables {u,v,t) 1-^ {—u,v,—t). Now if (0,05/02) is a saddle, then the 
singularities of X on §^ can be nodes or foci. Note that if 02—4(07 + 0507) > then 
{/± = 0}n§^, where f±{x, y, z) = y + 2o7z/ (— 02 + ■\/o2 — 4(07 + 05O7)), are invari- 
ant circles of X on §^ with cofactors K± (x, y, z) = (—02 ± \J a^ — 4(07 + a507))x/2. 
Therefore, in the case that the singularities of X on S"'^ are either saddles or nodes, 
then X has at least an invariant circle on Hence, by Corollary [21] X has no 
limit cycles on §-^. Thus, in this case, the phase portrait of X on the Poincare disc 
is equivalent to the one of Figures O (o) or (6) on §^ . 

Now, we consider the planar vector field induced by X through the stereographic 
projection (|10p at the point (0,0,1), i.e. we consider system ([Ti)) which in this case 
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Figure 5. Phase portrait of Case 3: ai = 04 = 0, 02 ^ and 07 ^ —05. 



becomes 

ii = — 2a5W + 4a2f ^ — 2(a5 + 2a7)u^f + 2a5t;'^, 
V — — 2a7U — 4a2uu — 2(07 + 2a5)uw^ + 2a7M'^. 

This system has an inverse integrating factor (see Remark II 7p given by V{u,v) ~ 
[u^ + v'^ + l)V2{u,v), where V2{u,v) = a7M'* + 2(a7(w^ - I) ~ a2v)u'^ + a'j{l + v'^f + 
2v[2a^v + a2(l — v"^)). Note that que the roots of V2 = in u are 



J a^^ajv^ + a2V + 07 ± K'|-\/"2 ^ 4(a^ + a^aj)) 

u = ±— . 

a? 

Therefore, if a\ — 4(a7 + 0507) < 0, V2 vanish only on the points (±1,0). Hence, 
by Remark [T7l system ((3T|) has a first integral defined on \ {(±1,0)}. Hence, 
system (|3T|) does not have limit cycles. Therefore, if the singularities of X on 
are foci, then the phase portrait of X on the Poincare disc is equivalent to Figure 
[51(c). 

Now if 07 — —05, then the system associated to X becomes 
i = 2/(a2y + 05^;), y = -a;(a2y + flgz), i = 0. 

Therefore, the phase portrait of X on §^ in this case is equivalent to Figure [4] (a) . 

Case 4- ai — a2 — 0. This case correspond to system ([25)) . We claim that we can 
reduce to the case 04 = 0. Now we prove the claim. In fact, consider the system 
associated to X on i.e. system ([3]) with ai = 02 = 03 = ag = 0, ag = —04 
and a\ + 0705 < 0. If 04 7^ and 05 + 07 7^ 0, by the orthogonal linear change of 
variables 

/a;\ ( a h Q \ ( X \ 

['i)-[t 

where a — a/i\/2a / {{a^ + a'j)-/6l3), b = {a5{a5 + aj)S— {2al + a5{a5 + ai))y/(T)/ ^/SPj, 
d = ~y/j 3/V2 S, e = -{ai{{a5 + a7)(4a| + 07(05 + 07) + 05(05 + 07)) + (07 - 
a5)V^))VWy, <^ = (05 + 07)2 + 4o|, cr = (05 + ajfS, (3 = 5- 7 = (05 + 
07)^(05 + 07 + 2o|) + (07 — 05) v^, we obtain 04 = 0. Here en denote the coefficients 
of system associated to X on the variables (a;,y, z). Note that in this case <5 > 0, 
cr > and as {5 - \/ct)(5 + ^) = 4o|5, ((05 + 07)^(0^ + o^ + 2o|) + (o^ - o^)^) • 
((05 + 07)^(01 + 07 + 204) - (07 - a|)v^) = 4(05 + a7)''(o4 + 0507)^, it follows 
that /3 > 0, 7 > 0. Moreover, we have that oi = 02 = 03 = 04 = og = og = 0, 
05 = — V2(o5 + a7)(o507 + a\)/^ and 07 = — ^/(V2(o5 + 07)). Note that 
0507 < 0. 
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Figure 6. Phase portrait of Case 4: ai = a2 = and 04 = 0. 



Now, if 04 7^ and aj — —05, then a| — 05 < and doing the orthogonal hnear 
change of variables 




04 + 05 



4 ^2(a5 +04)2 
04 + 05 



v/2(a5 + 04)2 

1 / 




we have as = 04(04 - aD/y/aKa^ + 04)^, 07 = \/al(a5 + 04)^/04, ai = 02 = 03 = 
04 = 06 = 08 = 0. Note that 050,7 < 0. In short, the claim is proved. 

Now we study the case 04 = and 05 + 07 7^ 0, X has six singularities on given 
by (±1,0,0), (0, ±1,0) and (0,0, ±1). The eigenvalues associated to the singular- 
ities (±1,0,0) and (0, ±1,0) are {0,±A/-(a^ + 05O7)} and {0, ±a/-(o^ + 05O7)}, 
respectively. As 0507 < 0, from ([30|) in this case if two of these opposite singu- 
larities are saddles the other two have complexes eigenvalues and they are foci or 
centers. In fact, they are centers. Suppose that — (05 ± 0507) < 0. By the change 
of variables (x, y, z) ^ (x, z, y) the system associated to X becomes 



(32) 



05ZJ/, ij = -(05 + 07)zx, i = ajyx. 



By Theorem [21 the singularities (0,0, ±1) are centers of the above system. Hence, 
(0, ±1, 0) are centers of X. Similarly if — (o^ + 05O7) < 0, by the change of variables 
(x, y, z) i—f (z, y, x), we conclude that (±1, 0, 0) are centers of X. Hence, the phase 
portrait of X on the Poincare disc is equivalent to Figure [6] 

If 04 = and 07 = — 05, then the system associated to X becomes 



X = osyz, y ■ 



0. 



Therefore, the phase portrait of X on in this case is equivalent to Figured (6). 

We note that the study of system ((32)) is called by some authors (see [6 ) the 
Euler Problem. 

Now we prove statement (d). If X has a semi-hyperbolic singularity then, by 
definition, for determining its phase portrait we have to study system ()15|) with 
04 + og 7^ and 0403 — 05O7 = 0. We can divided the proof in three cases. First 
case 04 = 05 = and so og 7^ 0. Second case 04 = 07 = and so og 7^ 0. Third case 
og = 05O7/04 and 04 7^ — og. However, we need study only the first case because 
the last two cases are equivalent to it by a rotation of the variables. In fact, if the 
coefficients of X satisfy the second or the third cases, then by the orthogonal linear 
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change of variables [x, y, z) — Mi(x, y, z), i — 1, 2, where 

f 

Ml = 



and 



V 



f- 



V 



as 




as 








as 







ai 







ay 




04 




04 


a| 


07 


a| 







V«7 + 




a| 





1 / 

o\ 



1 / 



the coefficients of X in the variables (5;, z) satisfy the first case, respectively. 
Thus, we suppose that 04 = as = and so as ^ 0. We distinguish the following 
four cases. 

Case 5: aiaj 7^ 0. We have that X has only four singularities on i.e. (0, 0, ±1), 
{ai/ yja^ + a\, 0, ai/ yja^ + a\) and {—a-j/ yja^ + af,0, — oi/i/a^ + a\). These sin- 
gularities correspond to the singularities (0,0) and (—07/01,0) of system ([15]). Do- 
ing the change of variables u = —asr/ar, v = r + s and introducing the new time 
variable r through rfr = —a%dt system ()15p can be written as 

(33) f = P2{r,s), s = s + Q2ir,s). 

If 0207 — oiog 7^ 0, we have that ip{r) = ((0207 — aia8)(a7-|-o|))r^/(o7a|)-f • • • is the 
solution of s -1- Q2{r, s) — 0. Hence, ip{r) = P{r, (f{r)) = (0207 — oias)r^/o§ + • • • . 
Thus, by Theorem 65 of page 340 of [1], we have that (0,0) is a saddle-node of 
system ([551) . Now, the trace and the determinant of the linear part from system 
([Is]) at (—07/01,0) are (0207 — ai08)/ai and 07(07 -|-of)/af, respectively. Therefore, 
(—07/01, 0) can be a node or a focus of system ([TS]) . Now, by the orthogonal linear 
change of variables 

/ ai 




1 

- 





07 





Ol 



the system associated to X becomes 




= P{x,y,z)^ 

= R{S:,y,z) = 



0207 — oios 

—xz - 



01O2 + a-ja% 

/ 2 1^*^^ 

y o^ + Oy y O]^ -r U7 

O1O2 +O7OS .2 , ~~ . 

— X -I- oixy 07x2, 

0207 - 



aryz, 



ails ~2 
-a; . 



V"? + «7 



Note that in these coordinates the singularity (— a7/(-\/af + af), 0, — oi/(-\/o^ + a^)) 
corresponds to the point (0, 0, —1). Hence, using the notation of Section[5l we have 
that R(u, v) — R{2u, 2v, + v'^ ~1) = 4(0207 — oiOs)/ + OyU^. Hence, R does 
not change sign and by Theorem [3l X does not have limit cycles on §^ . So the 
phase portrait of X on the Poincare disc is equivalent to one of Figure [Tj Note that 
the separatrices of the hyperbolic sectors of the diametrally opposite saddles-nodes 
can be connect or not, see Figures 5 and 6 of [Tl] . 
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Figure 7. Phase portrait of Case 5: aia-j ^ and 0702 — agOi 7^ 0. 




Figure 8. Phase portrait of Case 6: ay = and oi 7^ 0. 



Now if 0702 — agfli = 0, system p3p has a straight hne of singularities passing 
though the origin. Note that H(x,y,z) — ajx + aiz is a first integral of X and 
{H 0} n is a circle of singularities. So the phase portrait of X is equivalent to 
Figure 11 (6). 

Case 6: ay = and ai 7^ 0. So (0,0, ±1) are the unique singularities of X on 
Introducing the new time variable r through dr = —a^dt system (jlSp can be 
written in the form 

u = P2{u,v), v = v + Q2{u,v). 

We have that v — f{u) = — aiu^/2a8 + • • • is the solution of w + Q2{u,v) = 0. 
Hence, ip{u) — P{u, ip{u)) — a\v? /a\ + • • • . Thus, by Theorem 65 of page 340 of 
[T], we have that (0,0) is a topological node of system (fT5)) . 

Now we consider the planar vector field induced by X through the stereographic 
projection of (fTO]) at the point (0,0, 1), i.e. we consider system (fT4|) which in this 
case becomes 

(34) f — Aairs + 4a2S — Aa^rs , s = —2ass — Aair — Aa2rs + 208^ s — 2a8S . 

Since s \s=o— — 4air^, the orbits of vector field determined by system ([M]) intersect 
the straight line s — Q always in a unique direction except at the point (r, s) ~ (0, 0), 
i.e s = is transversal to the flow of system (|34l) except at the origin. Thus, 
X — {P,Q,R) does not have limit cycles on We also can apply Theorem [3] 
of Section [5] to give another proof of this fact. Then, it is sufficient to see that 
R{u, v) = R{2u, 2v, + v"^ — 1) = —Aagv'^ does not change sign. Hence, the phase 
portrait of X on the Poincare disc is equivalent to Figure [H 

Case 7: — and ai = 0. We have that {y = 0}n is a circle of singularities of 
X. As H{x, y, z) = asx + a2Z is a first integral of X , it follows that phase portrait 
of X is equivalent to Figure H (a) . 

Case 8: ay 7^ and oi — 0. Assume that 02 7^ 0, then X has only four singularities 
on S^, (0, 0, ±1) and (±1, 0, 0). The singularities (0, 0, ±1) correspond to the singu- 
larity (0, 0) of system ()15|) . Doing the change of variables u = —asr/a-j, v = r + s 
and introducing the new time variable r through dr — —agdt system (|15|) can be 
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Figure 9. Phase portrait of Case 8: 0702 7^ and oi = 0. 



written 



(35) 



r = P2{r,s), s ^ s + Q2{r,s). 



Since 02 ^ we have that s — (p{r) = (02(07 + af))r'^ / (araf) + ■ ■ ■ is the solution 
of s + Q2{r,s) — 0. Hence, ■>p{r) — P{r,Lp{r)) ~ (a7a2)r-^/a§ + • • • . Thus, by 
Theorem 65 of page 340 of [1], we have that (0,0) is a saddle-node of system 
([35l) . Now, the singularities (±1,0,0) have eigenvalues {0,— (a2± \J a^ — ^a^)/'2, 
and {0, (02 ± y^2|^^4af)/2, respectively. Therefore, (±1, 0, 0) can be nodes or foci 
oi X — {P, Q, R) restricted to Now, using the notation of Section O and doing 
the change of variables {x, y, z) 1 — > (z, y, x), we have that X becomes {V, Q, TZ) = 
{R{z, y, x), Q{z, y, x),P{z, y, x)) and so TZ{u, v) = P{v? + — 1, 2v, 2u) = 4a2W^. 
Hence, TZ does not change sign and by Theorem O X does not have limit cycles on 
So the phase portrait of X on the Poincare disc is equivalent to one of Figures 
ini Note that the separatrices of the hyperbolic sectors of the diametrally opposite 
saddles-nodes can be connect or not, see Figures 5 and 6 of pT] , 

Now if a2 = 0, X has a circle of singularities on §^ determined by the intersection 
of the plane a^x + a%y = with Note that H[x, y, z) = x \s & first integral of 
X. So the phase portrait of X is equivalent to Figure [4] (a). □ 

Proof of Corollary It is a straightforward consequence of the proof of Theorem 



7. Singular points of homogeneous polynomial vector fields of CP^ 

In this section we present arguments and results which go back to Darboux [3], 
see also [3]. First we recall that CP^ = {C^ \ {(0,0,0)}}/ ~ with the equivalence 
relation [X,Y,Z] ~ [X',Y',Z'] if, and only if, there exists A e C \ {0} such that 
[X',Y',Z']^\[X,Y,Z]. 

Let P, Q and R be homogeneous polynomials of degree m + 1 in the complex 
variables X, Y and Z . We say that the homogeneous 1-form w = PdX+QdY+RdZ 
of degree m + 1 is projective if XP + YQ + ZR = 0. Therefore, the 1-form is 
projective if and only if there exist three homogeneous polynomials L, M and N of 
degree m such that P = ZM-YN, Q = XN-ZL, R = YL-XM, or equivalently 



m 



□ 
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(P, Q, R) = {L, M, N) A {X, Y, Z). Then we can write 

L M N 
oj = X Y Z 
dX dY dZ 

We notice that the polynomials L, M and N are not uniquely determined by P, Q 
and R. These polynomials can be replaced by L' = L + XA, M' = M + YA and 
N' = N + ZA, where A is any homogeneous polynomial at the variables X, Y, Z 
with degree m — I. 

The projective 1-form uj defines a differential equation in CP^ given hy uj — 0, 
which may written as {ZM - YN)dX + {XN - LZ)dY + {LY - MX)dZ = 0. 

Usually in the literature a projective 1-form lu is called a Pfaff algebraic form of 
degree to + 1 of CP^ , see for more details Jouanolou JTT . 

Let F{X, Y, Z) be a homogeneous polynomial of degree n in the variables X, Y, 
Z and coefhcients in C. Then F{X, Y, Z) ~ is an algebraic curve ofCP^ Let 
p — (Xq, Yq, Zq) be a point of CP^. Since the three coordinates of p cannot be zero, 
without loss of generality we can assume that p = (0,0,-1). Then suppose that 
the expression of F{X, Y, Z) restricted to Z = — 1 is 

F{X, Y,~\) = Fd{X, Y) + Frf+i {X,Y) + --- + F„(X, Y), 

where < c? < n and Fj{X,Y) denotes a homogeneous polynomial of degree j in 
the variables X and Y for j — d, . . . n, with F^ different from the zero polynomial. 
We say that d — dp{F) is the multiplicity of the curve F = at the point p. If 
d = then the point p does not belong to the curve F — 0. If c? = 1 we say that p 
is a simple point for the curve F — 0. If d > 1 we say that p is a multiple point of 
_F = 0. In particular p is a multiple point of _F = if and only if 

dF dF dF 

Suppose that for d > we have that Fd = Ii\^iLY are different straight lines, 
called tangent lines to _F = at the point p, ri is the multiplicity of the tangent line 
Li at p. For d > 1 we say that p is an ordinary multiple point if the multiplicity of 
all tangents at p is 1, otherwise we say that p is a non-ordinary multiple point. 

Let F = and G = be two algebraic curves and p a point of CP^. We say 
that F — Q and G — intersect strictly at p, if F = and G = do not have any 
common component which pass through p. We say that F = Q and G = intersect 
transversally at p if p is a simple point of F = and G = 0, and the tangent to 
F = at p is distinct to the tangent to G = at p. The proof of the following two 
theorems can be found in [7|. 

Theorem 26. (Intersection Number Theorem) There exists a unique multi- 
plicity or intersection number lij), F C] G) defined for all algebraic curves F — 
and G — and for all point p of CP satisfying the following properties. 

(a) I{p, F n G) is a non-negative integer for all F , G and p when F ~ and 
G = intersect strictly at p. I{p, F (1 G) ^ oo if F ~ and G = do not 
intersect strictly at p. 

(b) /(p, FnG) = if and only if p is not a common point to F = and G = 0. 
/(p, F n G) depends only on the components of F = and G — which 
pass through p. 

(c) If T is a change of coordinates and T(j)) = q, then I(q^ T{F) D T(G)) = 

/(p,FnG). 
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(d) /{p,i^nG) = /(p,GnF). 

(e) I{p, FnG) > dp{F)dp(G) , the equality holds if and only if F ~ and G ~ 
do not have common tangents at p. 

(f) i{p, F n G1G2) = /(p, F n Gi) + i{p, F n G2). 

(g) I{p,FnG) = I{p, Fn{G+AF)) for all homogeneous polynomial A{X,Y, Z). 

Theorem 27. (Bezout Theorem) Let F = and G ~ be two algebraic 
curves of CP of degrees r and s respectively without common components. Then 
j:^I{p,FnG) = rs. 

Let i<i = 0, for i = 1, . . . , s, be algebraic curves and p a point of CP^. We define 
the multiplicity or number of intersection of the curves Fi ^ 0, . . . , Fs = at p as 
Hp, f^LlF^) = min,<,{/(p, n F,)}. 

A point [Xq, Yq, Zq] belonging to CP^ is a singular point of a projective 1-form 
Lu = PdX + QdY + RdZ of degree m + 1 if P{Xq,Yq,Zq) = Q{Xq,Yq,Zq) = 
R{Xq, Yq, Zq) = 0, i.e. [Xq, Yq, Zq] Satisfies the system 

ZM -YN = 0, XN - ZL^O, YL- XM = 0. 

The proof of the next two results can be found in j4j or [17] . 

Lemma 28. (Darboux Lemma) Let A, A' , B, B' , C and G' be homogeneous 
polynomials in the variables X , Y and Z of degrees I, V , m, m' , n and n' respec- 
tively, verifying the identity AA' + BB' + GG' = 0. Assume that the curves A — 0, 
B ~ 0, G = and the curves A' = 0, B' = 0, G' ~ do not have any com- 
mon component, respectively. Then I{p, A n B (1 G) + J2p Hp, A' f] B' D G') = 
(Imn + I'm'n')/^, where "f^l + l'^m-\-m'^n + n'. 

Proposition 29. (Darboux Proposition) For any projective 1-form of degree 
m+1 o/CP^ having finitely many singular points we have that its number of singular 
points taking into account their multiplicities or numbers of intersection satisfies 
J^p lip, [ZM - YN) n {XN - ZL) n {YL - XM)) + m + 

Proof of Theorem^ We have that X induce a projective 1-form uj on CP^ of degree 
n given hy u; = PdX + QdY + RdZ. Therefore, by the Darboux Proposition, w 
has (n— l)^ + (n— 1) + 1 = — n + 1 singular points, taking into account their 
multiplicities or numbers of intersection. Now, if [Xq, Yq, Zq] is a singular point of u; 
and {Xq, Yq, Zq) £ R^, then X{Xq, Yq,Zq), A S R, is a straight line of singular points 
of X which determine exactly two singularities of X on Hence, we can conclude 
that if X has finitely many singularities on S^, then X has at most 2(ri^ — n + 1) 
singularities on □ 

8. Topological classification of all phase portraits of homogeneous 
polynomial vector fields on §^ of degree 2 modulo limit cycles 

Let X — {P, Q,R) he a homogeneous polynomial vector fields on §^ of degree 2. 
Consider the system ([3]) associated to X. Using the notation of previous section, 
we can write ([3]) as 

X = z{a4X -\- a^y -\- a^z) — y{—aix — a2y) — zM — yN, 
(36) y = x{—aix — a2y) — z{—ajx — a^y — oqz) — xN — zL, 

z = y{~-ajx — asy — oqz) — x{a4X + a^y + a^z) = yL — xM, 
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where L{x, y, z) = ~ajx — a^y — a^z, M{x, y, z) = a^x + a^y + a^z, N{x, y, z) = 
~aix — 02?/, i.e. 

(37) 

Note that X = {P,Q,R) = iL,M,N)A {x,y,z). Then, (a;o,yo,^o) 7^ (0,0,0) is a 
singular point oiX if and only if there exist Aq € M such that (L, M, N){xo,yo, zq) = 
^a{xo,yo, zq). Now, as L, M and N are homogeneous polynomials, it follows that 
A(xo, j/o, Zq), X £ R, is a straight line of singularities of X. Therefore, for determin- 
ing the singularities of X we must find the real eigenvectors of the 3x3 matrix 
A given by ([57)) . i.e. we have to calculate the real eigenvalues of A. Let Aq S M 
an eigenvalue of A. If there exists a unique eigenvector associated to Ao, then it 
determines a straight line of singularities of X, and so we have two singularities 
of X on Now if there exist two eigenvectors associated to Aq, then we have a 
plane of singularities of X that determine a circle of singularities of X on S^. Thus, 
X has finitely many singularities on if and only if each real eigenvalue of A has 
exactly one eigenvector associated. Therefore, for degree 2, we obtain another proof 
of Theorem [T) We also conclude that if X has infinitely many singularities on , 
then it has one invariant circle on §^ formed by singularities and its phase portrait 
is equivalent to one of Figures H] or [T] 

In Section [6] we classified all phase portraits of homogeneous polynomial vector 
fields on of degree 2 having a degenerate singularity or a center. Therefore, 
for determining all phase portraits of a homogeneous polynomial vector field X on 
§^ of degree 2 we have to consider only the X's with all their singularities non- 
degenerate and without centers. In this case, the singularities can be saddles, nodes 
or foci. Now, as these singularities have index —1 or 1 and the maximum number of 
singularities of X is 6, by symmetry of X with respect to origin and the Poincare- 
Hopf Index Theorem (see, for instance [14]) it follows that we have to study the 
phase portraits with two saddles and four singularities that can be nodes or foci, 
and the phase portraits with two singularities that can be either nodes or foci. 

Let p be a singularity of X, by Proposition[T4l we can suppose that p = (0, 0, —1) 
and so as = ae — in Hence, 



A = 



We have that the eigenvalues of A are {0, (05 — ayi -y/a)/2}, where a = (05+07)^ — 
4a4a8. 

In the case that X has only two singularities that can be either nodes or foci, 
then A has either a unique eigenvalue with a unique eigenvector associated, or 
A has one real eigenvalue and two complex eigenvalue. In the first case is the 
unique eigenvalue of A and — (05 — aj + y/a){a5 — ay — y^) = 4(0408 — 0507), i.e. 
0408 — 0507 — 0. Hence, using the notation of Section [6l we have that (0,0, —1) is 
a degenerate singularity, because det{DX{0, 0)) = 0408 — 0507 = 0. This case was 
already studied in Section [H 

If A has one real eigenvalue and two complex, then is the real eigenvalue of 
A and a < 0. Hence, as in this case the system associated to X has the normal 
form ([3|) with = qq — 0, we have that R{u,v) — R{2u,2v,u^ + — 1) = 
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— 4(a4U^ + (05 + aj)uv + agw^). Solving R = with respect to u, it follows that 
u — —(05 + 07 ± y/a)v / {2ai) are its roots. Note that, as a < 0, 04 ^ 0. Thus 
we have that (0, 0) is the unique real root of i? = 0, and that R does not change 
of sign. Therefore, applying Theorem [3] (a), we have that in this case X has no 
limit cycles and so its phase portrait in the Poincare disc is equivalent to the one 
of Figure [H 

We have that (0,0,-1) is a saddle of X if and only if det{DX{0,0)) < 0, i.e. 
0408 — 0507 < 0. In this case, as (05 + 07)^ — Aa^as = (05 — 07)^ — 4(0403 — 0507) > 

0, A has three reals eigenvalues given by {0, (05 — 07 ± y/a)/2}. So a must be 
positive. Therefore, the eigenvectors of A associated to the eigenvalues {0, (05 — 
07± ^)/2} are (0, 0, 1), ((03(05 - 07 + ^))/(2oia8 - 02(05 + 07 + V")), (04(05 - 
a7 + ^))/(2o2a4-oi(o5 +07 - ^/a)), 1) and ((08(05 - 07 - ya))/(2ai08 - 02(05 + 
07 — ^/a)), (04(05 — 07 — ^/a)) / {2a2a4 — 01(05 + 07 + Va)), 1), respectively. Note 
that the eigenvectors of A are well defined if (2oia8 — 02(05 + 07 + ^/a)){2aias — 
02(05 + 07 - ^))(2o2a4 - 01(05 + 07 + ^))(2o204 - 01(05 + 07 - ^)) 7^ 0, i.e. 
if 160408(0304 — 0102(05 + 07) + 0108)^ 7^ 0. System l|36p with 03 = og = and 
04O8 — 0507 < have two saddles and four singularities that can be nodes, foci 
or centers if and only if A has three linearly independent eigenvectors. Thus, we 
distinguish the following cases. 

Case 1: 160403(0204 — 0102(05 + 07) + 0^03)^ ^ 0. This case has been classified in 
the previous paragraph. 

Case 2: 04 = 0, oi 7^ and 02(05 + 07) — 0103 7^ 0. In this case — 05O7 < 
and the eigenvalues of A are 0, — 07 and 05 with respective eigenvectors (0,0, 1), 
(0503/(02(05 + 07) - 0103), -05(05 + 07)7(02(05 + 07) - 0103), 1) and (07/ai, 0, 1). 

Case 3: 04 = oi = and 02 7^ 0. In this case —0507 < and so 05 + 07 7^ 0. The 
eigenvalues of A are 0, — 07 and 05 with respective eigenvectors (0,0,1), (1,0,0) 
and (0503/ (02(05 + 07)), -05/02, 1). 

Case 4- — ai — a2 — 0. In this case —a^aj < and so 05 + 07 7^ 0. The 
eigenvalues of A are 0, —07 and 05 with respective eigenvectors (0,0,1), (1,0,0) 
and (-03/(05 + 07), 1, 0). 

Case 5: 04 = 0, oi 7^ and 03 = 02(07 + 05)/oi. Now the eigenvalues of A are 0, 
—07 and 05 with respective eigenvectors (0,0,1), (07/01,0,1) and (—02/01,1,0). 

Case 6: oi = 02 = and 0403 7^ 0. The eigenvalues of A are 0, (05 + 07 ± ^/oi)/2 
with respective eigenvectors (0, 0, 1), (—203/(05 + 07 + ^/a), 1, 0) and (—203/ (05 + 
07 - ^/a), 1, 0). 

Case 7: 03 = 02(— 02O4 + 01(05 +07))/o^ and 01O4 7^ 0. The eigenvalues of A are 0, 
(02O4— oi07)/oi and (01O5— 0204)/oi with respective eigenvectors (0, 0, 1), (((01O7— 
02O4) (01(05 + 07) — a204))/(of (01(05 +07) — 202O4)), ((02O4 — 0107)04)/ (01(01(05 + 
07) — 202O4)), 1) and (— 02/01, 1, 0). Note that 04O8 — 05O7 — (0204 — 0i05)(ai07 — 
0204)/oi < and 01(05 + 07) — 20204 7^ 0. Otherwise, 05 = (202O4 — oia7)/oi and 
so 04O8 — 05O7 = (0407 — 0204)^/01 > 0. 

Case 8: 03 — 0. We have —a^aj < 0. This case is equivalent to one of Cases 2-5, 

1. e. if the coefficients of system ([55]) with 03 = og = satisfies 03 = then by the 
orthogonal linear change of variables 
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we obtain one of Cases 2-5. 

Case 1 is equivalent to Case 2 doing the orthogonal linear change of variables 
(a;, y, z) — M{x, y, z), where 



M 



( 



V 



2a8 



as + a? + \/a 



v/(a5 + 07 + V")^ +4a| V(«5 + 07 + ^Jaf + 4a| 
05 + a7 + -y/a 2a8 







^(05 + 07 + + 4a| ^(05 + 07 + ^JaY + 4a| 











1 / 



and the a^'s denote the coefficients of system p6p in the variables (x, y, z) satisfying 
Case 1. Indeed, denoting by hi the coefficients of system in the variables 
{x^y, 2), we have that ai = (2aia8 — 02(05 + 07 + ^Ja))/ \/ (05 + 07 + Va)^ + 4a|, 
02 = -(20208+01(05 + 07 + ^))/ (05 + 07 + ^JaY + 408, ^5 = 1/2(07-05 + ^"), 
07 — 1/2(05 — 07 + -^^), 08 = 04 + 08, and 03 = 04 = oe = 0. Now, as 160408(0204 — 
0102(05 + 07) + Oi08)^ 7^ 0, it follows that 04 = 0, oi 7^ 0, 02(05 + 07) — 0408 7^ 
and —0507 — a^a^ — a^aj < 0. 

Case 6 is equivalent to Case 4 if 05 + 07 7^ doing the orthogonal linear change 
of variables 

X \ /a b \ / X 
y ] = [ d e ] [ y 
z ) \ Q I ) \ z 

where o = ((a5+07)^(o4+08) + (o4-08)^/CT)/((o5+07)V2(^>/3), b = ( (05+07) (05(05 + 
aiY - 04(08(805 + 07) + 04(07 - 05))) - (04(04 - 08) + 05(05 + a^))^/a)/{^JZ|3^, 
d = ~^/]3/^/26, e = ((05 + 07)(o4(2o8(a4 - 08) - 07(05 + 07)) + 0508(05 + 07)) - 



(0704 + 0508)\/ct)/%/(5/37, 5 = (05 + a-j f + (04 - 08)^ 



(05 + 07)^0;, [3 



(05 + 07)^ + 204(04 — os) — -y/a, 7 = (05 + 07)^(0^ + 07 — 204O8) + (07 — o^^/ct, and 
the Oi's denote the coefficients of system (pS)) in the variables {x,y,z) satisfying 
Case 6. Note that in this case S > 0, a > 0. Now, as ((05 + 07)^ + 204(04 — 08) — 
y/a) ■ ((05 + 07)2 + 204(04 - 08) + V^) = 4a|(5, ((05 + 07)^(0! + 07 - 20408) + 
(of - al)^/a) ■ ((05 + ar fial + of - 20408) - (of - al)y/a) = 4(05 + a7)'^(o408 - 
0507)^, (05 + 07)^ + 204(04 — 08) = Og + of + 2a| — 2(0408 — 0507) > and 
Og + 07 — 20408 = (05 — 07)^ — 2(0408 — 0507) > it follows that f3 > and 7 > 0. 
Denoting by di the coefficients of system (I36|) in the variables {x,y,z), we have 

that 05 ^ -\/2(05 + 07)(0507 - 0408)/^, 07 = -y^/ (2(o5 + 07)), 08 = O4 + Os, 

oi = 02 = 04 = 03 = 04 = oe = 0. Moreover, —0507 < 0. 

Now, Case 6 is equivalent to Case 4 if 07 = —05 doing the orthogonal linear 
change of variables 



/ - 04O8 04\/o| - 04O8 + 05 



X 

y I = 

z 



0403 



O4 — 08 



y/aj- O4O8 



04 



V 







-y/ (05(04 - 08) + ^-(04 - 08)^0408)2 

"8 v^aj " a408 - 05 Vaj ~ "4fl8 

(05(04 - 08) + a/- (04 - 08)^0408)2 





X 

V 
z 



1 / 



Indeed, denoting by di the coefficients of system (1361) in the variables {x,y,z), we 
have that 05 = (05 + a804)(a8 — 04)/ (05(04 — as) + (04 — 08)^0408)^, 07 — 
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\J (05(04 — as) + ^J—{ai — ag)"^ a^as)^ / {a4^ — as), as = 04 + as and ai = 02 = 03 = 

0,4 = ae = 0. Now, since + a8a4 < 0, it follows that —0507 < 0. 

Case 3 is equivalent to Case 5 doing the orthogonal linear change of variables 



X 

y 



( 



\ 



as + ay 



yj(a^ +07)2 + a^ 
as + 07 

y(a7+&7)M~af 




\J (as + 07)2 + , 



V(a5 







1 / 




where the a^'s denote the coefficients of system p6p in the variables {x,y,z) sat- 
isfying Case 3. Indeed, denoting by hi the coefficients of system ((36)) in the variables 
(i, y, z), we have that ai = -a2(a5-|-a7)/ ^ (as + a7)2 + a|, 02 = -a2a8/V("5 + a?)^ 
5.5 = a7, 07 = as, as = 02(05 + aj)/di and 0,3 = 0,4 = de = 0. Now, as a2 7^ and 
— asa7 < 0, it follows that di ^ and —5,507 < 0. 

Case 7 is equivalent to Case 5 doing the orthogonal linear change of variables 
(x, y, z) = M (i, y, z) with M given by 



V 



0204 — 01(05 + 07) 



0401 



a/ (oi (as + 07) - 0204)2 + a^a 
0401 

y/ (oi (os + 07) — 0204)2 + a^o 




^y (01(05 + 07) - 0204)2 + a2a| 
0204 - 01(05 + 07) 

^(01(05 + 07) - 0204)2 + a^ol 










1 / 



where the a.^'s denote the coefficients of system ([5^ in the variables {x, y, z) satisfy- 

(0204 — aia5)(oia7 — 0204)701 < 



-05O7 

-,2^2 



ing Case 3. Note that, as oi ^ and 0408 - 

0, it follows that (01(05 + 07) — 0204)2 + a2a4 7^ 0. Denoting by di the co- 
efficients of system p6p in the variables (a;,y, z), we have that oi = (01(05 -I- 
07) — 20204)01/ (01(05 -I- 07) — 0204)2 -I- a2a|, 02 = (04(02 — Oi) — 0102(05 -I- 
"7))/ V (ai(a5 + 07) - 0204)2 -I- a\al, 05 = (0204-0405)701, 07 = (0204-0107)701, 
Os — 02(05 -I- 07)701, and 03 = 04 = og = 0. Now, as 04(05 + 07) — 202O4 7^ 0, it 
follows that oi 7^ 0. Moreover, —0507 — (0204 — oia5)(oi07 — 0204)701 < 0. 

We conclude that for determining all phase portraits of all homogeneous poly- 
nomial vector field on §2 of degree 2 having two saddles and four singularities that 
can be nodes, foci or centers we have to consider only Case 2, 4 and 5, i.e. we have 
to consider the following families 

X = 05J/2 - ?/(-oix - 02?/), 
y = a;(-aia; - 02y) - z(-07a; - osy), 
z = yi-ajx - asy) - a^xy, 

' aj) — oiOs 7^ and —a^aj < 0; 



(38) 

with oi 7^ 0, 02(05 



(39) x = b^yz, y ^ -zi-bjx ~ bsy), z ^ yi-bjx - bgy) - b^xy, 

with —6567 < 0; and 



(40) 



X = c^yz - y{-~cix - C2y), 

I X / C2(C5 + C7) 

y = x{~cix - C2y) - z \-c^x y 

( C2(C5 + C7) \ 

z = y\-c-jx — yj-c^xy, 
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(a) {b) (c) 



Figure 10. Homoclinic orbits. The symbol • denotes a focus, a 
center or a node 

with Ci 7^ and — C5C7 < 0. 

Consider the vector field X determined by system (|39|) . Doing the orthogonal 
linear change of variables 



Note that for 02 ~ ^bs, 05 = — (65 + ^7) and 07 = 67 system (j4T|) becomes the 
system of Case 3 in the proof of Theorem [S] of Section [S] with + a^a^ < 0. 

Now we will study the possible connections of separatrices of families ([55)1 and 
(|40| on We have consider only the homoclinic orbits, because as the saddles 
of families ((38|l and (j40| on are diametrally opposite, by Proposition \T2\ if X 
has a heteroclinic orbit on then X has an invariant circle on Then, the 
phase portraits of homogeneous polynomial vector fields on of degree 2 having 
invariant circles have been classified in [IV . 

We saw that families (p8)) and (|40| have six singularities, two saddles and other 
four singularities that can be nodes, foci or centers. Therefore, the homoclinic 
orbits that can appear in families (l38|) and (|40|) on §^ must be one of the described 
in Figures [10] (a) and (c) . 

Proposition 30. Let X be a homogeneous polynomial vector field of degree 2. Then 
X cannot realize Figure (b) or (c) . 

Proof. Suppose that X realizes some of Figures [TU] (&) or (c). We can suppose 
without loss of generality that the three singularities of Figure [10] (6) or (c) are in 
the southern hemisphere. We consider a great circle on §^ that contains the saddle 
and another great circle that contains the other two singularities. The great circle 
through the saddle is chosen in a such way that the two straight lines determined 
by these great circles through the central projection ^ in the tangent plane to §^ 
at the point (0,0, —1) are parallel. Eventually these two great circles can coincide. 
Now, there exists a different great circle C such that the straight line determined 
by it is also parallel to the previous straight lines, it does not contain the saddle 
point and the other two singularities belong to the half-sphere determined by it 
and containing the saddle point. As the orbits of X are symmetric with respect 
to the origin, these three invariant circles cannot intersect in two singular opposite 
points. Moreover, we can choose C such that the straight line determined by it is 
sufficiently close to the straight line that contains the saddle point, and C intersects 
transversally the homoclinic orbits in four points (see Figure [TT|) . Note that, as the 
orbits of X are symmetric with respect to the origin, C has six tangent points with 




system (|39)) becomes 
(41) X = -bsy- 



„-;2 



{h + b7)yz, 



y = bsxy + b7xz, z = b^xy. 
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orbits of X. But this is a contradiction with Proposition [231 Therefore, the phase 
portrait of X cannot reahze Figure [10] (6) or (c). □ 

From Figure [HI we can see that there exist homogeneous polynomial vector fields 
on §^ of degree 2 having a homoclinic orbit of the type described in Figure [TOj (a) , 
where the singularity inside the loop is a center. However, we do not have examples 
where this singularity is a focus or a node, but numerical computation show that 
they exist. For example, we can consider the following family of homogeneous 
polynomial vector fields on of degree 2 determined by the system 

(42) X = —xy — 2y^ + yz, y = x^ + 2xy + xz + a%yz, z = —2xy — a^y'^ , 

with as 7^ 4. System (|42|) is obtained from ([SS]) substituting ai = —1, a-z = —2, 
a5 = 1 and ay = 1. Now, consider the system induced by (I42|l through the central 
projection on the tangent plane to §^ at (1,0,0), i.e. system ([9]) becomes 

, , ii — — u — — 2u^ — 2u^u — oguu^, 

■i) = — u — agf + + 2uu — 2uu^ — agw'^. 

This system has three singularities, (0,0), (1,0) and (^—asjias — 4),2/(a8 — 4)) 
with respective eigenvalues (— ag ± Y^a| + 4)/2, (2 — ag ± \J{2~ a^Y — 16)/2 and 
(ai - 2a8 + 2 ± ^{a^ - 2as + 2)2 - 8{al + (ag - 4)2 + 4))/(2(4 - ag)). Using the 
software PA, we can obtain the phase portraits of system in the Poincare disc. 
Thus, for ag = 8/5 and ag — 9/5 we obtain the phase portraits given by Figures 
[T2](a) and (6), respectively. In Figure [T^we have that (0,0) is a saddle and (1,0), 
(— ag/(ag — 4),2/(ag — 4)) are unstable foci. Note that in Figure [T^ (a) one of 
the stable separatrices of (0,0) emerges from the unstable focus (1,0). Now, in 
Figure [T^] (6) one of the unstable separatrices of (0, 0) tend to the unstable focus 
(1,0). Therefore, there exists at least one limit cycle surrounding (1,0). We can 
conclude by Figure [T^j that there exists a value of the parameter ag belonging to the 
interval (8/5,9/5) such that there is a homoclinic orbit surrounding the unstable 
focus (1,0). In this homoclinic orbit borns a limit cycle because in Figure [T2l (6) 
the point (1,0) is an unstable focus surrounded by an unstable separatrix coming 
from the saddle (0,0). It seems that moving the parameter of the system this limit 
cycle ends in the Hopf bifurcation which takes place at the singular point (1, 0). 
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Figure 13. Phase portrait on the Poincare disc of homogeneous 
polynomial vector fields on §^ of degree 2 with a saddle. The 
symbol • denotes either a focus or a node surrounded perhaps by 
some limit cycles or a center. Note that the three singularities also 
can be surrounded by some limit cycles. 

We can conclude that if X is a homogeneous polynomial vector field of degree 
2 on §^ with all its singularities non-degenerate, then its phase portrait on is 
topologically equivalent to one of the phase portraits on given in Figures 1131 
This proves Theorem [8l 

9. HOPF BIFURCATION ON HOMOGENEOUS POLYNOMIAL VECTOR FIELDS ON §^ 

OF DEGREE 2 

Let X = {P,Q, R) he a. homogeneous polynomial vector fields on §^ of degree 
2 and let p be a singularity of X. Without loss of generality we can suppose 
that p = (0, 0, —1), because we can do a rotation of S0{3) which preserves all the 
properties of X. Hence, using the notation of Section [51 the singularity p is a weak 
focus of X on §2 if ti{DX{Q, 0)) = 0, det(L»X(0, 0)) > and p is not a center. 

Proposition 31. Let X be a homogeneous polynomial vector vector field on §^ 
and let p = (0,0,-1) a singularity of X . Suppose that the system associated to X 
is in the form ^ with — oq = 0. Then p is a weak focus of X if and only if 
as ~ —cii, 04 + 0705 < and 04(0^ — 03) + 0102(05 + 07) ^ 0. Moreover, p is stable 
if 04(01 — 03) + 0402(05 + 07) < and unstable if 04(01 — 03) + 0102(05 + 07) > 0. 

Proof. The proof of this proposition follows directly of the proof of Theorem [2j □ 

Before stating the next theorem, we will state the Hopf Bifurcation Theorem. 

Theorem 32. (Hopf Bifurcation Theorem) Suppose that the analytical parame- 
trized system x = X{x, p), x G R^, /i € R, has a singular point at the origin for all 
values of the real parameter fi. Furthermore, suppose that the eigenvalues, Ai(/i) 
and X2{fi) of DX{0, fi), are pure imaginary for p = fj,Q. If the real part of the eigen- 
values, Re\i{fi) — ReX2{fi) in a neighborhood of fiQ, satisfies — (i?eAi(/i)) |^=^j,> 
and the origin is an asymptotically stable fixed point when ^ — fiQ, then 

(a) fi = /io is a bifurcation point of the system; 

(b) for p e (/^i,^o); some < /io, the origin is a stable focus; 

(c) for /i G (/io, ^J^2), some p2 > Mo; the origin is an unstable focus surrounded 
by a stable limit cycle, whose size increases with fi. 

The proof of Hopf Bifurcation Theorem can be found in [15' . 
Proof of Theorem\^ Consider the vector field X = {P,Q,R) of the beginning of 
this section. By Sections [HI and [51 it follows that if X has a weak focus, then the 
system associated to X is equivalent to one of systems or ([IO)) . 
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Consider system pOj) with ci 7^ and — C5C7 < 0. 
gularities, (0,0, ±1), {cr/{y/4+l?^),0,ci/{^/^T^)), 



This system has six sin- 
i-cr/{^ 



c 7),0,-c i/( 

\/cl + cf ) ) and two singularities on , given by (— C2 / ( a/c| + ) , ci / ( ^/c^'+l^ ) , 0) 
and (c2/(-\A^f+cf), — ci/(-\/cf+cf ), 0) with respective eigenvalues {0, (\/c| + cf ± 



4(c2 + C5C7))/2}, {0, (-y^fT^±y^ 



4(ci 



C5C7))/2}. ASC5C7 



> 0, the singularities of X on are either nodes or strong foci on S^. 

The planar system determined by system (|40p through the central projection on 
the tangent plane to §^ at the point (0, 0, —1) is given by 



2 / , N 2 C2(C5 + C7) 2 
U — —C5V + ClUV + C2V — (C5 + Ci)U V uv , 



-C7U - 



C2(C5 + C7) 



Cl 



Cl 



U — CiU^ — C2UV ~ (C5 + cj)uv'^ 



C2(C5 + C7) 3 



Cl 



see This system has two singularities (0,0) and (— C7/ci,0) corresponding to 
the singularities of X that does not belong to S^. We saw, in Section [51 that (0, 0) 
is a saddle. Now, (— C2C5 ± \/ c^cl — 4(07 + C5C7)(cf + c'y))/2ci are the eigenvalues 
of (— C7/ci,0). Again, as C5C7 > 0, a necessary condition in order that (— C7/ci,0) 
be a weak focus is C2 — 0. If C2 — 0, by the orthogonal linear change of variables 




system (j40|) becomes 

X = -ciy"^ + cjijz, y ^ cixy + C5XZ, z = -{c^ + cr)xy. 



This system satisfies Case 3 in the proof of Theoreni[5] Therefore, {—'^il \J'A + ^7, 0, 
— ci/(-\/c^ + C7) is a center of X on §^ and the phase portrait of X on is equiv- 
alent to the phase portrait given by Figure [5] (&) or (c). Thus, in the family ([^0]) 
its singularities cannot be weak foci and so it does not have a Hopf bifurcation. 

Now we consider system (|38|) . System (pS)) has six singularities, (0,0, ±1), 
[a^ai|^fo,-a^(a^\■a'J)|^fa, {a2{a' i^aT)- a\a%) j ^ fo), {-a ^as/^/^, 05(05 +07)/^^, 
-(02(0 5 + 07) - QiQ8)/\/g), [aTl^a^i \- a^,0, oi/a/o^ -|- of ) and {-a-jj^a^ \- af,0, 
—axj^Ja^ -|- of), where cr = 05(0^ -f (05 -1-07)^) -I- (02(05 -1-07) — oiog)^. We saw in 
Section [8] that (0,0, ±1) are saddles and the other four singularities can be nodes, 
foci or centers. By the orthogonal linear change of variables 






Ol 



system ((38)) becomes 



X = — 



\Ja\ + o2 
07 

\ yof+o? 



0102 -I- 0708 _ _ 





07 



Ol 



X 

y 



0207 



OiOg^^ _2 

==^a;z - oiy 



aiyz, 



ai02 -f 0703 .2 , ~~ , f , \ — 
— a; aixy + (05 -f 07)a;z, 



0207 - 



OiOs ^2 
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Note that in these coordinates the singularity {—aj/{y/af+(^),0, —a\l{\J a\ + a^)) 
goes over the point (0,0,-1). Let ax = — (aia2 + a-ja-g) / \/a\ + Oy, 5,2 = — oi, 



04 = —(0207 — aiag, / ^ a\ + , 05 = —07, 07 = 05 + 07, and as — 0. By 
Theorem [21 {—a-j/^Joj + af,0, — ai/i/a^ + a^) is a center of system ([55]) on 
if a4 = 0, i.e. 0207 — aias = 0, 04 + 0507 = —(07 + 0507) < and a4^(a\ — 
a|) + 01(12(05 + 07) = -\/afT"afai05ag/o7 = 0. However, this is not possible, be- 
cause 010507 ^ 0, 0207 — oios = (\/af + afai0508)/a7 ~ implies 02 = os = 0, 
which is in contradiction with the hypothesis 02(05 + 07) — oiog ^ 0. Therefore, 
{—O'l I \/ + a?: O7 — oi/ -y/oy + Oi) cannot be a center. Now, if 0207 — oiOs = 0, by 
Proposition [31] it is a weak focus of system p8)) . 

The planar system determined by system (|38|) through the central projection on 
the tangent plane to at the point (0, 0, —1) is 

, > il — —051; + OiUU + 021"^ — (05 + 07)1*^11 — OgUU^, 

(44) . 2 / \ 2 3 

V = —a'ju — ogw — oiu — 02UW — (05 + a-j)uv — ogu , 



see This system has three singularities (0,0), (—07/01,0) and 

(—0508/(02(05 + 07) — oiOs), 05(05 + 07)/(o2(o5 + 07) — oiog)) corresponding to 
smg ularities of X on S^. The singularity (—07/01,0) corresponds to singularity 
(— 07/y^a7 + af, 0, — oi/y^of + of) of system (|38)) on §^ and its eigenvalues are 
'^±(/^) = (Ai± _ 4(a2 ^ a^a-!){a\ + of))/(2oi), where ^ = -(0207 - oiOg). We 
saw that if ^ = and \/ cl\+ 0^010508/07 < 0, (—07/01,0) is a stable weak fo- 
cus. Moreover, for — 2-\/ (07 -I- 0507)(a^ + o^) < ^ < 2-^/ (07 -|- 0507)(a^ + o^) the 

singular point (—07/01,0) is a strong focus and — — (ReA±(/i)) |n=o= 1- Hence, by 

O/i 

the Hopf Bifurcation Theorem, we have a Hopf bifurcation. This proves that there 
exist examples of homogeneous polynomial vector fields on §^ of degree 2 with at 
least one limit cycle. 

Now, consider the orthogonal linear change of variables {x^y,z) — M{x,y,z), 
where M is the matrix 

/ 05-1-07 08(02(05 -I- 07) — O1O8) 0508 \ 

v/(o5 -1-07)2 +al ^(05 -1-07)2 +aly/a 

0.8 (05 -I- 07)(02(05 -I- 07) - OiOg) 05(05-1-07) 



^/{a5 +aj)^ +al a/(o5 + aj)^ + ajy 



la 



\ 







05 vTas + OTp+of 02(05 -I- 07) - 0108 



Hence, system ([38)) becomes 



la 



y = —ciiSp' — a2Xy + drxz, 
z = — a42;2 — (5g + a^)xy, 

where oi = — (oi02((o5 + a'jY — ai) + CLs,{a^a1 + (05 -I- a7)(a2 — a\ + 05(05 -I- 
07))))/ (V(a5 + 07)2 +o|^), 02 (02(05 -I- 07) - 0i02)/v/(a5 + 07)2 -|-a|, 04 = 
-(08(0207 - O8O1) - 01(05 + a^ay))/^, 05 = 05, and 07 = -(05 + 07). Note 
that in these coordinates the singularity (—0508/-ya-, 05(05 -t- 07 )/\/ct, —(02(05 -I- 
07) — aitts) / y/a) becomes the point (0,0,-1). By Theorem [21 (— 0508/y^, 05(05 -I- 
07) / \/a, — {02(0,5 + 07) — Oi08)/Vct) is a center of system (1551) on §2 jf = 
0, a| + 0507 < and 04(0,1 — Oj) -|- 0102(05 + 07) = 0. We have that 04 — 
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implies that 08(0207 — ogOi) — ai(o| + 0507) = 0. Note that in this case 
os 7^ 0, otherwise 04 = — ai(a| + a^a7)/^/a 7^ 0. Hence, 02 — 01(05 + 0507 + 
a|)/(o7a8) and a\ + 0,50.7 — — (05 + 0507) < 0, 04(01 — 03) + 0102(05 + 07) = 
aiyJal{ol + (05 + 07)2)(a2(a5 + 07)2 + ol{o\ + o^))/|a7a8|. As oi ^ 0, (-osog/ 
y/a, 05(05 +07) — (02(05 +07) — oiogj/v^) cannot be a center. Now if 08(0207 — 
a80i) — 01(05 + 0507) = 0, by Proposition [3T] it is a weak focus of system ([55)) . 

The singularity (-0503/(02(05 + 07) - ai08), 05(05 +07)7(02(05 + 07) - oiog)) of 
system (|44l) corresponds to the singularity {—a^a^/^/a, 05(05 + oj)/ y/a, —{02(05 + 
07) — 0i08)/v^) of system ([55]) on §^ and its eigenvalues are A±(^) — (fi ± 
i//L<2 - 4(05 + a507)cr)/(o2(o5 + a7)-aia8), where /i = -(a8(o207-080i)-oi(o| + 
05O7)). We saw that if ^ = and oi/|o708|-\/(a^(o5 + 07)^ + o|(o^ + Oy)) • 
y^ajial + (05 + 07)2) < 0, then (-0503/(02(05 + 07) - 01O8) , 05(05 +07)/ (02(05 + 
07) — O1O8)) is a stable weak focus. Moreover, for — 2-\/ (05 + 0507)0- < /i < 
2V(«5 + 0507)0-, (-0503/ (02(05 + 07) - 0103), 05(05 + 07)/ (02(05 + 07) - 0103)) 

is a strong focus and — (ReA±(/i)) L=o= 1- Hence, by the Hopf Bifurcation Thco- 

rem, we have a Hopf bifurcation. Note that system p8|) cannot have two weak foci 
at the same time. □ 



10. Rotated vector field family 

In this section, we summarize the behavior of the limit cycles in the special one- 
parameter family given by a rotated family of planar vector fields. The earliest 
work about these families can be found in the paper [5] of Duff in 1953. Later on 
Seifert T8j, Perko [Hj and Chen Xiang-yan [19, 20, 21], etc successively improved 
the work of Duff. 

Consider the vector fields Xa{x,y) = {P{x,y,a),Q{x,y,a)) depending on the 
parameter a. Suppose that when a varies on an interval (o, 6), the singular points 
of the vector fields Xa remain unchanged, and for any fixed point p = {x,y) and 
any parameters ai < a2 (z (o, 6), we have 



(45) 



P{x,y, ai) Q{x,y, ai) 
P{x,y,a2) Q{x,y, a2) 



> (or < 0), 



where the equality cannot hold on an entire periodic orbit of Xa with a = ai, for 
i = 1, 2. Then, the family of vector fields X^ is called a / generalized) rotated family 
with respect to the parameter a. Here, the interval (o, b) can be either bounded or 
unbounded. 

The geometric meaning of condition ()45p is the following. At any fixed point 
p — (a;, J/), the oriented area between the vectors {P{x,y,ai),Q{x,y,ai)) and 
(P(a;, y, a2), Q{x, y, a2)) has the same (or opposite) sign as sgn(Q;2 — ai). That is, at 
any point p = (a;, y), as the parameter a increases, the vector (P(a;, y, a), Q{x, y, a)) 
can only rotate in one direction; moreover, the angle of rotation cannot exceed tt. 

In the following we present four important results concerning periodic orbits and 
limit cycles for rotated vector field families Xa- 

(i) Non-intersection property: For distinct ai and q;2, the periodic orbits of 
the vector field Xa with a = ai and of the vector field Xa with a = a2 
cannot intersect. 
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(ii) Stable and unstable property: When the parameter a changes shghtly and 
monotonically, the stable and unstable limit cycle cannot disappear; it ex- 
pands or contracts monotonically. 

(iii) Semiestable property: When the parameter a varies in the suitable direc- 
tion, a semistable limit cycle bifurcates into one stable and one unstable 
limit cycle. When a varies in the opposite direction the semistable limit 
cycle disappears. 

(iv) Ending or starting property: When the parameter a varies a family of limit 
cycles only can disappears or appears either at a singular point, or in a 
semistable limit cycle, or in a separatrix cycle, or at infinity (i.e. the family 
becomes unbounded). 

Consider system (|^n|) with ci 7^ and — C5C7 < 0. This system has six singu- 
larities, (0,0, ±1), (C7/ Vcj + C^, 0, Ci/y^cj + C^), (-Cy/y /cj + c'^, 0, - Ci/y^cj + C^) 

and two singularities on S^, given by {—cil\/ci^ + , ci/ \/c\ + c\,Q) and (02/ 
■y/c2 + c\, — Ci/-\/c| + c\,Q). We saw in the proof of Theorem [H] that the singu- 
larities of X on §^ are either nodes or strong foci on 

Now we consider the planar system determined by system (j40p through the 
central projection on the tangent plane to §^ at the point (0, 0,-1), i.e. 
(46) 

2 ( , \ 2 C2(c5 +c^) 2 
u — —c^v + Cluv + C2V — (C5 + C7)u V uv , 

Cl 

C2(C5+C7) 2 / , X 2 C2(C5+C7) 3 

V = —ClU V — ClU — C2UV — [C5 + Cj juv V , 

Cl Cl 

see dni). This system has two singularities (0,0) and (— C7/ci,0) corresponding to 
the singularities of X that does not belong to S^. We saw, in Section [51 that (0, 0) 
is a saddle and, by the proof of Theorem [9l if C2 = then (— C7/ci,0) is a center; 
and if C2 0, then it is a node or a strong focus. 



Proposition 33. The vector field associated to system (j46[) does not have limit 
cycles. 

Proof. Consider the one parameter family of vector fields (P(u, w, C2), Q(w, w, C2)) 
associated to system (|46l) . We have that 

,~ ^ C5(C5 + C7) 2 , 2x 2 

(C2 - C2) (1 + U +V )V . 

Cl 

Therefore, system (l46l) determines a rotated vector field family with respect to the 
parameter C2. Now, note that for C2 = 0, (— C7/ci,0) is a center of this family. 
Then, by the non-intersection property (i) and the ending or starting property 
(iv), system (^S)) does not have limit cycles. □ 



P{U,V,C2) Q{U,V,C2) 
P{U,V,C2) Q{U,V,C2) 



We conclude by Proposition [33] that if the vector field X on §^ associated to 
system (i40l) with ci ^ and — C5C7 < has a limit cycle, then this limit cycle 
surrounds one of the singularities that belongs to §^ . These singularities always are 
nodes or strong foci. We believe that the following conjecture holds. 

Conjecture The vector field associated to system ([iO)) with Ci ^ and — C5C7 < 
does not have limit cycles on 
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